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Nomenclature 



Vector flux per unit length 
Magnetic flux density 
Electric flux density 
Electric field intensity 
!!agnetic field intensity 
Electric current density 
Magnetization 
Area 

Permittivity of free space 
Dielectric constant or permittivity 
Electric charge density 
Conductivity 

Permeability of free space 
Relative permeability 
Resistivity 


A 

B tesla * w/m^ 

0 

D coulomb/m^ , 

E volt/m ; 

H amp ere /m * 

3 ampere./m^ \ 

ft ampere/m j 

A m^ t 

\ 

e Q farad/m \ 

^ 1 

i 

z 1 

< 

p a c/ m3 j 

i 

o mhos/m 

y 0 henry /m j 

i 

y 

p ohm-m 


\ 



! 
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I. Introduction 


Cryogenic turoomachinery of the type used to pump high pressure 
fuel (liquid H 2 ) and oxidizer (liquid C 2 ) to the main engines of the 
Space Shuttle have experienced rotor instabilities. Subsynchronous 
whirl, an extremely destructive instability, has caused bearing failures 
and severe rubs in the seals (1,2). These failures have resulted in 
premature engine shutdowns or, in many instances, have limited the 
power level to which the turbopumps could beoperated. The labyrinth 
seals originally used in these pumps were initially indicated as a source 
of subsynchronous vibration (2). Other principal sources of self- 
excited instabilities in the hydrogen pump, in addition to the seals, are 
aerodynamic cross coupling turbine and impeller forces and internal 
shaft hysteresis and friction forces caused by relative motion between 
surfaces. All of these mechanisms can induce self-excited rotor non- 
synchronous whirl motion in a pump (3-6). The hydrogen pump, for 
example, has all thre.e instability mechanisms present because of its 
built-up strucrure of spline fits and high energy density level (7). 
The SSME oxygen pump as well as the hydrogen pump is susceptible to 
self-excited whirl motion (8). 

The occurrence cf self excited instability can be extremely dan- 
gerous because the whirl amplitude of motion may increase rapidly with 
increasing energy input. Unlike synchronous vibrations whose ampli- 
tudes reduce as the critical speed is traversed, subsynchronous whirl 
orbits may spiral out until metal to metal contact occurs in the impellers 
and seals. The occurence of metallic rubs on the oxygen pump is 
particularly serious as catastrophic fires may occur (9). High syn- 
chronous vibrations in a pump may be controlled by proper balancing 
and by avoiding operation near the critical speeds. However, with 
self-exciteo' whirl motion, improvement of balance has little or no effect. 
In fact, it may even aggravate the situation (10). 
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The large roior orbits caused by self-excited whirl induce high 
be/ ing loads. Since rolling element bearing life varies approximately 
.iversely. as the third power of loading, an increase in hearing loading 
can lead to a dramatic reduction in bearing life. Current turbopump 
designs do not include provisions for multiplane trim balancing of the 
built-up rotor after final assembly in the pump casing. Although trie 
impellers, turbine wheels, seal runners and the shaft may be individ- 
ually balanced, a satisfactorily balanced assembly is not always guar- 
anteed. 

The need for dissipating vibrational energy in high performance 
turbomachines has long been recognized. Many of today's turbojet 
engines in both civilian and military aircraft incorporate vibration 
dampers at or near bearing supports (11). 

With the availability of bearing lubricating oils in aircraft engines, 
a device known as a "squeeze film" damper has been used quite success- 
fully in attenuating potentially large and damaging forces (11). Over 
-the years, numerous investigators have produced both analyses and test 
results on squeeze film dampers (12-21). Because of these efforts, it 
is possible to design such dampers with the assurance that they will 
perform reliably in man/ different applications. Viscous "squeeze film" 
dampers have successfully attenuated both synchronous and nonsyn- 
chronous whirling if properly designed. 

In the Space Shuttle Main Engine (SSME) turbopumps, liquid h 2 
and LOX are used to cool the rolling element bearings. Because of the 
extremely low viscosity of the liquids (liquid H 2 has a viscosity approxi- 
mately equal to air at room temperature), they cannot be considered as 
adequate in providing a damping media fcr either viscous ^hear or 
squeeze film damping. Unless suitable energy dissipating devices can 
be developed, future generation turbopumps may be susceptible to the 
same potentially destructive vibrations as have been encountered in the 
current generation of cryogenic turbomachinery (9). 

The objective is to examine one of the damping mechanisms that might 
be suitable for the development of a practical discrete cryogenic machinery 
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damper. Listed below are the more common damping mechanisms available. 


1. Viscous shear and squeeze film bearings 

2. Visco-elastic material dampers, such as rubber isolation pads 

3. Coulomb-friction dampers 

4. Turbulent flow close clearance seals 

5. Eddy-current or magnetic dampers 

The first and probably most common damping mechanism is viscous 
damping; here the damping force is directly proportional to velocity. 
The constants of proportionality differ however, in the case of "squeeze- 
film" damping, as compared to viscous shear damping. "Squeeze-film" 
coefficients vary directly with viscosity and as the cube of damper 
length while varying inversely as the cube of the clearance. To obtain 
any effective damping with very low viscosity fluids such as liquid H 2 
and 0 2 either the size of the damper must be made quite large and/or 
the clearance between moving and stationary members should be made 
quite small. The damping coefficient for a shear film damper is directly 

proportional to the fluid viscosity and the shear area and again in- 

versely proportional to the first power of the clearance. This type of 
damper would provide very little damping unless made prohibitively 
large and would not be acceptable in a compact turbomachine. 

A second type of damping, visco-elastic or hysteresis damping, is 
generally produced by elastic materials. Almost all materials exhibit 
some sort of damping when strained repeatedly. Visco-elastic materials 
such as rubber, and used in machinery, provide a sizeable amount of 
damping and are generally quite effective. Obviously rubber would not 
be suitable at cryogenic temperatures since it would lose its visco- 

elastic properties, i.e., become very brittle (22).. 

A third type of damping used in machinery isolators is friction or 
coulomb damping. In a friction damper the force is directly propor- 
tional to the coefficient of friction of the contacting surfaces, the area, 
and the pressure applied to bring the plates into contact. The damping 
force in this case is not proportional to the velocity. The problem 
with using this type of damper in rotating machinery is the inability to 
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predict the amount of damping available for any given situation. Values 
of the coefficient of friction are unreliable. The contacting surfaces 
under too little pressure, slip relative to one another or, with too great 
a pressure, do not move at all and therefore povide little, if any, 
damping. A considerable effort was expended by Rocketdyne to incor- 
porate a colomb-friction damper into the SSME hydrogen pump bearing 
supports. This effort was unsuccessful. 

The stability of the hydrogen pump was eventually improved by 
incorporating close clearance seals and stiff bearing supports based on 
design recommendations of the Rotor Dynamics Laboratory of the Uni- 
versity of Virginia. The turbulent flow seals produce both principal 
and cross-coupling stiffness and damping coefficients. Under proper 
selection of bearing support stiffnesses and seal clearances, the seal 
effects can promote rotor stability. However, critical speeds are now 
placed in the operating speed range, and when clearance seal wear 
occurs, this stabilizing effect is lost. 

A fifth damping mechanism is the eddy-current or magnetic damper. 
Many devices based on this type of damping are currently being used. 
Most of these applications are in instruments where the damping forces 
required are quite small. 

The damping force Is velocity or frequency dependent but more impor- 
tantly, the damping coefficient varies Inversely as the resistivity of a 
conductor moving in the magnetic field. If such a damper were to be used 
in a liquid Hg pump, for example, the extremely low temperatures encount- 
ered would significantly decrease the resistivity of the conductor 
material, thereby producing a reasonable value of damping. 

This report outlines the efforts of a preliminary study of the feasi- 
bility of using an eddy-current type of damping mechanism for the SSME. 
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II . Fundamental Design Equations for Eddy Current Damper 


The fundamental electromagnetic field equations , 
pertaining to the design and analysis of an eddy current 
cryogenic pump damper, are derived from Maxwell's equa- 
tions . 

These equations are as follows: 


2.1. Electro-Static Field Equation . The displacement 
current density J io related to the charge density o 2 put 
unit volume by 


J 



» 


where $ is the average velocity. 

The electric current across a surface S is defined 
as the rate at which charge crosses that surface. 

The current flow across this surface is given by 


f s ■ [ s - 


ndS 


The 

the charge 


gradient of the displacement field is equal 
density 


to 


- p 8 

where D is the density of electric flux passing through 
a given area. This relationship is the differential 
form of Gauss' Law. 


( 2 . 1 . 1 ) 


( 2 . 1 .?) 


(2.:. 3) 


2 . 2 Electro-Magnetic Field Equation - Ampere's Law 
(Maxwell's First Law) . The ampere- turn drop around a 
closed circuit equals the current enclosed 

1-9 H‘3l 
* -x 

where tt - magnetization vector. 
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From Equation (1.1.2), 

1 ’ { 5 ' dI ' f[ h 


dS 


By Stokes Law 


£ H*dl - JJ" (V x H)-dS 


Hence V x H 


s 

-* 

J 


(2.2.2) 


(2.2.3) 


By replacing 3 by the total current density which 
is the stab of the conduction current density and the 
displacement current density, then Equation (2.2.3) can 
be written as 

7 X 8 - £ + f! . (2-2.4) 

The above is called Maxwell's first equation. 

2 . 3 Faraday's Law of Induction (Maxwellb Second Law) . The 
Faraday law of induction states that the voltage V 
(induced) is equal to 

V induced - -N 

Where N ■ number of turns in coil 
♦ flux linked 

Let ^ be called the electric field and be defined 
as the gradient of a voltage (scaler function) by 

$ - W . 


(2.3.1) 



(2.3.2) 
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The voltage is defined by the contour integral 
of the electric field by 


V - ^ $.3l . 


By S toko's theorem 

V = <£ E-dl = ff (v x E).dS 

«■ JJ S 

From Equation 1.3.1 for one turn 

V * ‘3t * 

If we define the quantity $ as the flux density 
or magnetic induction, then 


$ 


V • 


We can write 


$ (v X E)dS - -f- £ ff ».dS . 

S JJ S 

Dropping the integral sign, we have 

'«* ■ -If • 

This is the Maxwell's second equation. 

2. 4 Divergence of the Magnetic Field (Maxv - 111b Third 
Law) . If we consider the flux passing through a pie- 
shaped section radial to a current I in a conductor, 
then the flux entering the section also emerges from the 
section and the net flux build-up is zero. 


or 


* enter ~ ^emerge " ® 

B *dS • B . dS 

enter emerge 


0 , 


(2.3.3) 


(2.3.4) 


(2.3.5) 


(2.3.6) 


(2.3.7) 


(2.4.1) 
(2 4.2) 


i .t •: 


2-3 


© 


but Equation 2. 4. 2 represei.ts a closed surface integral, 
^ B*dS * 0, 

and by use of the divergence theorem , 

B*dS = Xjy* (v*B) du = 0 

or (v*B) = 0 . 

This is Maxwell's third equation and simply states that 
magnetic fields neither emerge from nor close at a point. 


/ ^ /. ON 

\A.t. J/ 


(2. 4. A) 
(2.4.5) 


2 . 5 Constitutive Relations . To the above differential 
equations are added the constitutive relationships des- 
cribing the macroscopic properties of the medium being 
dealt with in terms of permittivity s , permeability y 
and conductivity o. The quantities e, u, and o are not 
necessarily simple constants. For example, in the case 
of ferromagnetic materials, the B-H relationship may be 
highly nonlinear. These constitutive relations are given 
by 


I - 

el 

(2.5.1) 




B - 

UH 

(2.5.2) 

-*■ 

-► 


J - 

oE . 

(2.5.3) 
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2 . 6 S implified Equations for an Eddy-Current Damper * 
Magnetic Induction - E? 


The flux (j>g in webers for a magnetic field can be defined in exact 
analogy with the flux for the electric field, namely 

$6 = J ^ / (2.6.1) 

where B is the basic magnetic field vector called the magnetic induction 
in gauss or webers/meter 2 and the integral is taken ver the surface 
for which is defined. 

The definition of f is as follows: If a positive test charge q Q is 

fired with velocity v through a point P and if a (sideways) force £ acts 
on the moving charge, a magnetic induction B is present at point P, 
where £ is the vector that satisfied the relation 

? = q Q v x 3 , (2.6.2) 

v, q Q , and £ being measured quantities. The magnitude of the mag- 
netic deflecting force ?, according to the rules for vector products, is 
given by 

F = q Q v B sin 0 (2.6.3) 

where 0 is the angle between v and if. The magnetic field always acts 
at right angles to toe direction of motion. 

The units of B are: 

i tesla = 1 weber/m 2 = IQ 4 gauss = 1 nev/ton/coul(m/sec) = 1 newton/amp-m 

Faraday's Law of Electromagnetic Induction 

Faraday's Law of Induction says that the induced emf 4 ' n a cir- 
cuit is equal to the negative rate at which the flux through the circuit 
is changing. If the rate of change of flux is in webers/sec, the emf 4 
will be in volts. In equation form, 

{ * - • ( 2 - 6 - 4 ) 

The minus ,'ign is an indication of the direction of the induced emf. 

Lenz's Law states that the induced current resulting from the induced 
emf will appear in such a direction that it opposes the change that 
produced it. 

*From material and relations found In D. Halliday and R. Resr.ir.k, Physics, 

John Wiley & Sons, Inc., 1965. 
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As an example, consider Figure 1 which shows a rectangular 
loop of wire of width L, one end of which is in a uniform magnetic field 
£ pointing at right angles to the plane of the loop. This fif.ld of & 
may be produced in the gap of a large permanent magnet or an electro- 
magnetic. The dashed lines show the assumed limits of the magnetic 
field. The experiment consists of pulling the loop to the right at a 
constant speed v. The flux <|>g enclosed by the loop is 


* B = BLx 


(2.6.5) 



Figure X. A rectangular loop is pulled out of a magnetic 
field with velocity v. 

From Faraday's Law the induced emf is 

£ . . = . SiSL*! * . Bl^S = BLv, (2.6.6) 

* at at oi 

where - 4rr was made equal to the speed v at which the loop is pulled 
dt . 

out of the magnetic field. This induced emf sets up a current in the 

loop, determined by the loop resistance R, 

i = i s . (2.6.7) 

1 R R 

From Lenz's Law, this current must be clockwise in the above figure 
since it opposes the change (the decrease in q>g) by setting up a field 
that is parallel to the external field within the loop. 
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ORIGIN"- ^AG£ IS 
OF POOR QUALITY 


The current in the ioop will cause forces F 2 , and F 3 to act on 
the three conductors, as given by equation 2.6.3. Because ? 2 and ? 3 
are equal .and opposite, they cancel each other. which Is the force 

that opposes any effort to move the loop, is given in magnitude from 
equations 2.6.3 and 2.6.7 as 

r2i 2w 

Fj = i L B sin 90° = ^ — (newton). 


For completeness, a check of equation 2.6.8 for units is as follows: 


i (amp) 


1 

R 


.■ volts -. _ BLv . weber 
^ohm ' R meter 2 


x meter x 


meter 1 
sec x ohrrr 


BLv . weber - 
R ' sec-ohm'' 


or 


ohm 


, weber * 
' amp-sec' 


Thus, 


e - B 2 L 2 v 
R 


. weber 2 
' meter 4 


x meter 2 


meter amp-sec * 

sec weber * 


- B 2 L 2 v weber-amp 
R ^ meter *' 


and since 


then 


1 (weber) = ( 


newton -meter 
amp 


), 


B 2 L 2 v . newton-meter amp . 
R ' amp * meter ' 


or 

F = B - 2 - L ~ (newton) 


The resistivity, p, is a characteristic of a material rather than of 
a particular specimen of a material and has the units of ohm-meter. 
Since the resistance of an electrical conductor is directly proportional to 
the length of the conductor and is inversely proportional to the cross 
sectional area, it is related to the resistivity as follows: 


( 2 . 6 . 8 ) 
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Thus, 

yields 


R = p ^ (ohm-meter *gggs) (2.6.9) 

substituting this relation for the resistance R into equation (2.6. 8) 


F = v (newton) 


( 2 . 6 . 10 ) 


This force, ?, is directly opposite to the direction of motion of the 
conductor and could be considered as the eddy-current damping force. 


It is noted that the eddy-current damping force is directly pro- 
portional to the velocity, F » v, and by introducing a constant of 
proportionality called the damping coefficient, C^, then 

? = C d v (2.6.11) 

It is obvious, then, that the damping coefficient would be 


r - B 2 LA .. newton-sec . 
^d p ' meter ' 


( 2 . 6 , 12 ) 


where B = magnetic flux density in webers/meter 2 , L ~ length of con- 
ductor in meters, A = cross section of conductor in meters 2 , and p * 
resistivity of conductor material in ohm-meters. 


The mks system and the English system of units for the damping 
coefficient are related as follows: 

i newton -sec _ . newton -sec lb meter in -3 lb-sec 

1 meter ~ meter 4.448 newton x 39.37 in ” ” inch 


1 QWgfflaK , 5.710 x 10' 3 

meter inch 


lb-sec . . in 2 newton -sec 

inch ' ’* 31 x ,u meter 
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Sample Calculation - Damping Coefficient 

Assume a circular 4 inch diameter copper conductor with a cross- 

2 

sectional area of 0.0625 Inch is immersed in pressurized liquid H 2 at 

27°K and the magnetic flux density is 7,C00 gauss or 7 x 10 1 webers/ 

-9 

meter 2 . A value of p = 2 x 10 ohm-meters for copper is very conser- 
vative. For a high purity copper at 27°K, the resistivity can be as low 

-9 

as 0.14 x 10 ohm-meters.* 




_ B 2 LA _ ( 7 X 10" 1 ) 2 x (n 10.2 x 10 ~ 2 ) x (0.64 10 ~ 2 ^ 


2 x 10 


(2.6.13) 


r = 3220 = 18 . 4 

d meter inch 


• ! 
\ » 


Sample Calculation - Stiffness (k) 

Assume the same constants and conditions as above and determine 
k when the peak-to-peak amplitude of motion, is 0.005 inches at a whirl 
frequency of 190 Hz. The stiffness k is 


_ af _ B 2 LA de _ B 2 LA ew , newtons 


ae 


pe dt 


pe 


C 1 


meter 


) 


(2.6.14) 


-i 2 
r\ * \ t 


,-2 


-2^ 


x-4 


- (7x10 'f x n(10.2x10 c ) x (0.64x10 ) 1.27x10^ Q . 

— n a X o ^ * I 


2 x 10" 9 x 1.27 x 10 " 4 


k * ’•» X 106 ■ 11 ' 000 i£h 


» : * # * . 
t * - . 


» . » , * 


* R. Barron, Cr yogenic Systems , McGraw-Hill Book Co., 1966. 
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t 
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III. FINITE ELEMENT GA1EKKJ.N ANALYSIS OF EDDY-CURRENT LOSSES USING 
TWO-DIMFNSlONAL POTENTIAL VECTOR FORMULATIONS 

INTRODUCTION 

In this report, the method of solution of the eddy-current problem 
is presented. The report presents two formulations used in the solutions 
of two-dimensional eddy-current problems: These formulations are called 

the magnetic vector potential method (MVP) and the electric vector 
potential method (EVF) . Both methods generate a second order Helmholtz 
type complex differential equation. The formulation presented here is 
taken from J. M. Schneider in his Ph.D. Thesis on "The Finite Element- 
Boundary Integral Hybrid Method and Its Application to Two-Dimensional 
Electromagnetic Field Problems", R.P.I., 1982, under the direction of 
Dr. Scheppard J. Salon, Department of Electrical Engineering, Rensselaer 
Polytechnic Institute, Troy, New York. 

For the case of the eddy-current damper moving in a magnetic field, 
the EVP formulation is applicable rather than the MVP formulation. The 
partial differential equation is multiplied by a weighting function and 
the Galerkin ijethod is applied to generate a finite element approximation. 
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M 


0 


3.1 Derivation of Maxwell's Sinusoidally Time-Varying Field Equations 


, t t t . 3D 

1. ? x R « J c + aF 

or in integral form 

/tvit .^(J c+ |2)m s 

The first Maxwell equation is referred to as Ampere's Law and 
H ■ magnetic field. 


(3.1.1) 


(3.1.2) 


2. VxE--ff 


or in integral form 


ft) 


The above equation is called Faraday's Law. 


(3.1.3) 


(3.1.4) 


3. 7-D ■ p 

or in integral form 

£ D*ds ■ J pdv (Gauss's Law) 

4. 7«B "0 or 6 B’ds ■ 0 
(nonexistance of monopole) 

Note that the point and integral forms of the first two 

equations are equivalent under stokes theorem, while the point 

and the Integral forms of the last two equations are equivalent 

under the divergence theorem. For free space, where there are 

no charges (p ■ 0) and no conduction currents (J £ « 0), Maxwell's 

point equations assume the following form (1) : 

* * 

* * A 
•> -► 

7*D “0 

7*B “ 0 


(3.1.5) 


(3.1.6) 


(3.1.7) 


(3.1.8) 


(3.1.9) 

(3.1.10) 

(3.1.11) 


0 


However, reduction of Equations 1.1 - 1.4 to the sinusoi- 



dally time-varying magnetoquasi-static equations requires the 
following assumptions (2) : 

1. All fields vary sinusoidally with time. 

2. Displacement currents and surface charges 
are neglected. 

3. Free charges and surface currents are non- 
existent . 

Assumption No. 1 was assumed by Schneider for power appara- 
tus which operate at relatively low frequencies. The assumption made 
that all field quantities vary sinusoidally from D.C. to several 
hundred hertz is appropriate, since it is the eddy-current pheno- 
mena in the sinusoidal steady state which is to be modeled. 

The second assumption on the displacement current where 
3D 3 ± 

3F e 3F E (3.1.12) 

is neglected is due to the relatively low permittivity (e) of most 
material, and can be neglected in the presence of the conduction 
current J, Hence, the magnetic field and conduction current are 
assumed to be predominant in the operation of an eddy-current 
damper * 

Maxwell's sinusoidally tims-varying magnetoquasi-static field 
equations are: 


V x E « 

-j“B 

(3.1.13) 


-*■ 


7 x H - 

J 

(3.1.14) 




V • B « 

0 

(3.1.151 

V • J * 

0 

(3.1.16) 


j 

1 

» 

i 
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In addition to Maxwell's Equations are the constitutive 


relations: 


B = pH 


J - aE 


(3.1.17) 

(3.1.18) 


and the interface conditions 


n • (B x - b 2 ) * 0 


(3.1.19) 


n x (H^ - H 2 ) - 0 

-v -+■ -*■ 

n x (E x - e 2 ) * 0 


(3.1.20) 


(3.1.21) 


The effective damping generated by eddy-currents in a 
conductor moving in a magnetic field is related to the power 


loss which is given by: 


. xsj lii 2 


v a 


(3.1.22) 


A quantity of use in the computation of inductance is the 
magnetic energy W^. For nonlinear materials it is expressed by: 

»■%/// B*Hdv 

m y 


(3.1.23) 


For linear-permeable materials. Equation 1.23 reduces to 


“» • *///„ 


(3.1.24) 


C 
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3.2 Magnetic Vector Potential Formulation 
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The Maxwell' 8 equations cannot be utilized to solve the problem 
of eddy-current losses in a conducting sheet moving in a magnetic 
field. The use of vector potential functions has been employed 
by numerous authors to reduce the coupled first order partial 
differential equations into a single second order Helmholtz wave 
equation. The basic two approaches which Schneider presents in 
detail are called the magnetic vector potential (MVP) method and 
the electric vector potential (EVP) method. 

The formulation of a two-dimensional eddy-current magnetic 
field problem requires the utilization of a vector potential func- 
tion, rather than a scalar potential function. The reason for this 
is the inability of a scalar potential to adequately describe the 
vector properties of the eddy-current generated in a conducting 
sheet . 

The magnetic vector potential method is based upon the use of 
a single component vector whose curl is the magnetic flux density 


given as follows: 


V x A - B 


V • A ■ 0 


Using Maxwell's equations and the constitutive relationships 


and the vector identities 


7* (7 x V) - 0 


7 x VU - 0 


where V and U are arbitrary vector and scalar functions, 


respectively. 


7 x(7 x E) - -7 2 E + 7(7*E) - -7 2 E 
if 7»E - 0 


© 
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The following partial differential equation for the magnetic 
vector potential is developed 



where J q is the applied current density in the Z direction, 
vector A is assumed acting in the Z direction. 

The expression for the resultant current density vector J is 
composed of the applied current J q and the eddy-current reaction 
juoA. 

J ■ J q -jiocrA 

The expressions for rhe magnetic flux density components in 
the X and Y directions are given by 

i j h 

3 _ 3 _ 3 _ 

3X 3Y 3Z 

0 0 Az 

B - — i - ~ A 


B - V x A - 


Hence, the magnetic flux density components are assumed to lie 
in the X— Y plane and are normal to the magnetic potential vector. 
The power loss is given by 

Mi dv .Jjf M 

neglecting the applied current J q and considering only the eddy- 
currents generated 

Je - jtooA 



(3.2.6) 


(3.2.7) 


(3.2.8) 


(3.2.9) 


(3.2.10) 


(3.2.11) 
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The eddy-current pjwer dissipated Is given by 


P 


€■ 


2 <J (A A*) dv 


(3.2.12) 


For a two-dimensional conductor with an effective skin 
penetration depth of 6e, the power loss is given by 


P 


(ir o 6e 


(A A*) dX dY 


(3.2.13) 


-'-'a 

where A* * complex conjugate magnetic potential vector function. 

The MVP approach considers only one current component and 
two V.vx components. This approach is used for the computation 
of the magnetic field distribution, inductance and eddy-current 
losses in any power apparatus having one predominant current 
component. Inherent in the single component MVP approach is 
the existance of coupling between current carrying regions insula- 


ted from one another. 



3.3 Electri c ..or Potential Formulation 

'.’he general eddy-current problem to ad > .. its most general 
form is extx difficult because it r--q». . . . j. three-dimensional 

analysis of the He.i.<tthor.::z aquation i ' Lon, if the conductor 

is saturated, then the problem is non as well. The problem 

of eddy-current generation is often reduced in complexity by con- 
sidering simplified cases where the structure is either long or 
planar. In the first case, the electric field and current density 
possess only one component. 

In the second case, the magnetic field is assumed to have 
only one component, while the current function may have two com- 
ponents. 

The magnetic vector potential formulation is used in the 
first cas^ while the electric * ector potential method is used for 
the planar representation. In both cases, the resulting partial 
diff e. dntial equations are similar in nature to the general 
Helmholtz equation. 

To employ the electric vector potential approach, Schneider 
assumes an electric potential function EVP (A) such that the curl 
of the EVP function is equal to the current density. 

Thus, 

V x A ■ J 

If A ■ A(z)k only, then 

7 * A ■ (if ) 1 - (if) j 

where 



(3.3.1) 


(3.3.2) 


(3.3.3) 


3-8 



From Maxwell's equation relating H and J 




0 ' 


* i 


i l 
f 

\ ; 
\ * 


f 

-* f- 

- t 

. i 




;i 


7 x K » J 

Hence, the magnetic field strength vector H must be different 
from the EVP vector A by an arbitrary vector. Since the curl of 
tl/.. vector must vanish, the arbitrary vector must be equivalent 
to the gradient of a scalar function. Thus, we have 

4 4 * 4 1 

H - A - V* (3.3.4) 

If Ho is the excitation magnetic field intensity, then 

H - Ho + A - V(j> (3.3.5) 

Following the procedure of Carpenter and selecting the 
-*■ 

Coulomb gauge for A such that 

-+■ -+■ 

V • A « 0 

eliminates the sources of $ and removes the V$ term from 
Equation 3.5 

H - Ho + A (3.3.6) 

Assuming that all field quantities are independent of the 
Z coordinate and that the permeability and conductivity are 
linear and Isotropic, the following partial differential equa- 
tion for A ™ A(Z) only is obtained 

33? + ayf ” “A ■ a^Ho (j.3.7) 

where 

o 2 ■ jwyo 

The resultant magnetic field intensity H is composed of an 
excitation component Ho and an eddy-current component: A with the 
corresponding current density components given by 


t 3A 
J X ’ 3Y 


j - _ M 

• Y 3X 


as stated in Equation 3.3. 


c 
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( 3 . 3 - 8 ) 


In Maxwell’s equations, across an interface, the condition 
n x - H 2 ) - 0 

must be met- 

This condition may be fulfilled by requiring the normal 
derivative of A to be discontinuous by the ratio of region 
conductivities 

3Ai 01 '3T 2 

3 n 02 an 

where T) is the unit vector normal to che interface. 

The seccnd or EVP approach allows for the generation of 
two-dimensional currents. In actuality, the eddy-current problem 
is three-dimensional, but the approximations of the two-dimensional 
solution is reasonable and is governed by the effective depth of 
penetration 5. 

Inherent in the single-component EVP approach, as stated by 
Schneider, is the absence of coupling between conductive regions 
insulated from each other. The reasons for this is that the 
induced eddy-currents require a conductive path in order to flew 
from one region to another and also that only the excitation 
agnatic field incident to a particular conductive region is 
modified by the eddy-currents induced into that region. For the 
case of the eddy-current damper analysis, the EVP "ormulatton will 
be utilized. 

The solution to an EVP(A) problem can be represented graphi- 
cally in the form of an equipotential plot. The difference of 
potential between two constant A lines equals the total per unit 
depth current flowing tangentially between them. Hence, A is 
a current describing function. 
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3 . 4 Galerkin Finite Element Formulation of the Electric Vector Potential 


Boundary Value Problem 

The vector form of the Helmholtz equation may be solved by a 
finite element proximation. In the majority of finite element 
formulations, a variational principle must first be obtained. 

The governing partial differential equation is transformed into 
an equivalent integral or functional statement of the problem. 
Minimization of the functional yields the desired solution. 
Finlayson and Striven state as early as 1967 that there is no 
practical need for variational formalism.. The Galerkin method 
or the method of weighted residuals is straight forward and avoids 
completely the effort and mathematical embellishment of a varia- 
tional formulation, and that apart from self-adjoint linear 
systems, which are comparatively rare, there is no practical need 
for variational formalism. 

Although the EVP partial differential equation is linear, 
self-adjoint and amenable to variational formalism, the Galerkin 
approximation is straight forward and is preferred. It is of 
interest to note that for any linear self-adjoint partial differ- 
ential equation, the variational and Galerkin finite element 
approximations yield identical simultaneous equations. 

Consider the general Helmholtz equation of the form 

3 ( 1 3a\ 3 (l 3a\ . A _ 

3X ^y ax ) 3Y (y 3 y) ” -*“ oA " F 

where 

A * Unknown potential function varying sinu- 
soidally in time with frequency. 

y,o ■ Material properties. 

F ■ Known excitation function. 
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Let H be an approximate solution to Equation 4.1 in a planar 
region* Since A will not in general satisfy Equation 4.1, a 
residual R will result. 


B _ 3 /l 3 a\ ^ d /l 3A \ . . — 

R 3X dXJ + 3Y 3Y J J A “ F 

In order to minimize the residual error R, it is multiplied 
by a weighting function W and the weighted Integral over the 
region is set equal to zero. 



0 


A major advantage of the Galerkln method is that it can be 
used to reduce the order of the derivations in the partial differ- 
ential equation. This will have a considerable significance on 
the choice of the shape functions as to the permis ible order of 
the polynomial functions used in the analysis. 

The Galerkin weighted integral Equation (4.3) is first 
transformed by the Divergence theorem and Green's theorems which 
are obtained from the Divergence theorem. 

For an arbitrary vector 


Iff 


V«Adv 


-y -*• 
A*nds 


In cartesian coordinates 


*-f £* 


+ A cos(y,n) + A cos(z,n)] ds 
y z 


cos(x,n) 



(3.4.2) 


( 3 . 4 . 3 ) 


(3.4.4) 


(3.4.5) 


(3.4.6) 
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Expanding the volume integral and rearranging 


uV 2 Vd\ 


r* * r 

- I Vu* VVdv +| u 
* v s 


BV 

3n 


dS 


The above equation is known as Green’s first identify. 
Green’s second identity is given by 


(UV 2 V-W 2 u) dv 



3V ,. 3u . 

“ m- v to. 1 oS 


dR 


The expanded form of Equation (4.3) Is 

+ h (vlf| 

- jo) WaA dR -fr WF dR « 0 

R — R 

Applying Green's first identity 


- I f lfwi£ . 3W 3a\ fwBA 

J I y I^X 3X + 3Y 3Y ) dR+ J y Sr. 


ds 


- J/« 




i' 


WaA dR - / / WF dK 


(3.4.7) 


(3.4.8) 


(3.4.9) 


(3.4.10) 


The contribution due to the contour integral is assumed zero and 
thereby implicitly '.atisrying the homogenous Newman boundary conditions. 
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The second order partial differential equation is reduced to the 
following integral equation with first order derivations 




3A 

ax 


+ 


aw 3A 

3Y 3Y 


dR 




WF dR 


In the Galerkin method of weighted residuals, a set of shape 
function? may be chosen for the whole domain. These functions must 
also satisfy the boundary conditions. 

In the first order finite element approximation of Equation (4.11), 
the region R is divided into a number of triangular elements over each 
of which A varies linearly and y and a are constant and the forcing 

A 

function F is taken as its average value F. 

Applying these assumptions. Equation (4.11) becomes 



The summation extends over the y elements into which region R 
ie divided. In the first order finite element method, each subregion 
m is represented by a triangle having local nodes 1, j, k,at its 
vertices. The values of the function A at the nodes 
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(3.4.11) 


(3.4.12) 


(3.4.13) 
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For any point within the triangular region, the value of A is 


<B 


given by 


A - CN ± Nj N k ] ( A, 


(3.4.14) 


The functions N^(X,Y) are called shape functions and have the 


property that 


N i (X i’ V " 1 


N^, ) - N 1 (X k , Y k ) - 0 


The value of A within the triangle is assumed to be linear of 


the form 


A(X, Y) - a + bx + cy 
m 


(3.4.15) 


Solving for A^ in terms of the coordinates Y^ of the point 


« a + b x^ + cy^ 
■ a + b Xj + cy^ 
- a + b x k + cy k 


m element 


(3.4.16) 


Solving for the shape functions N^, N ^ , N k we obtain 


(a^j, + bj x + c i y)/2Am 
(a^ 4- bj x + c^y)/2Am 
(®k + b k x + c k y)/2Am 


(3.4.17) 


fc; 

iV-- 

► * . 
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Where 
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a i 

“ X j \ - Y j \ 

S j 

“ *k T i- Y k X i 

a k 

“ X i Y j - Y i X j 

b i 

" Y j - Y k 

b J 

- y - y 
k i 

b k 

“ Y i- Y j 

C i 

“ 

C j 

m x i~h 

C k 

« Xj - x t 

Where Am * area 

of triangle, 

- <X 3 - X,) « k - V - (X k - 


■ b i c j - b y c i ’ \ C J- b j 
’ b k c l- b i c k 


The shape functions given by Equation (4.17) correspond to element m 
and only to points within the boundary region. The linear shape functions 
take on the value of unity at a node i and vary linearly to zero at the 
opposite nodes j and k. The value of at the node i for element m is 
the same as the corresponding point j for element N. Although the 
function A is continuous in the first order finite element method, the 
first derivative is not continuous in passing from one element to an 
adjacent node. 
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The Galerkin finite element approximation is developed by choosing 
the shape functions as the weighting functions for element m. 


W - 



(3.4.18) 


The partial derivatives are given by 




1 

24m 


b i b i 


1 

2 Am 





(3.4.19) 


(3.4.23) 


For the first order finite element approximation, the derivatives 

A and A are constant in the element m. 
x y 


aw _ 

3X " 2Am 



aw _ _i_ 
3Y " 2 Am 



C k 


(3.4.21) 


(3.4.22) 
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The first Integral 


i_ /IE 3A 3W SAjdR 

ym I 3X 3X &Y 3Y / 
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[ b * b ; <] 

K \ 


+ C 


[ c i b J c k] 


Since the matrix coefficients are constant, the integral reduces 
to Am. Expanding we obtain 


& 


b. J + c .2 


_A_ b i b j + c i c j 
Am4ym 

b i b k + c i c k 


b b + c c 
i j c i j 

b 2 + c 2 

j j 

b. b. + c. c, 
j k j k 


b b + c c 
i k i k 


Vk + <: j c k (3 -‘- 24 > 

\ ! + s 2 


The integral 


becomes 


juo / WA dR 
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jlOO 


N 



OF 

[2 

N N, 

i 

i j 

L N. 

N 2 

i 1 

i 

[, N 

N N, 

i k 

j k 


N i N k 
N J N k 

V 



A i 

dR 

A 3 

w J 


Carrying out the required integrations, we obtain 


j too Am 
12 


2 

1 

1 


1 

2 

1 


1 

1 

2 


— 

A. 


i 

# 

A. 


J 

J 

w ^ - 


The following function integral is given by 
dR 


- Fm 




-Fa 


X" 

dR ■ - ~ Am 

r i i 

N 3 

»— i 



LiJ 


3A 

Let ipm ■ ^ along a boundary element. 
The contour integral is given by 

Vtymdc ■ - — f Wdc 

pm 


■/" 


/ 


( 3 . 4 . 25 ) 


( 3 . 4 . 26 ) 


( 3 . 4 . 27 ) 


( 3 . 4 . 28 ) 


The normal derivative of A, i|<m is assumed to be constant over each 
segment on c and is independent of A. Let it be a coordinate along c^ 
with its origin at node j and directed towards node k. 


■V (x k ' V 2 + <Y k - V 2 


Let L 
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The single-component EVP pattlal differential equation for linear 
isotropic materials is 

|p + - ° 2 A - a 2 Ho - F (3.4.30) 

where 

a 2 ■ jtu pa 

Ho * applied magnetic field intensity 


The finite element equations for the m** 1 element is given by 


A 

S ii 

S 13 

S ik 


- ^ 

A 

l 


r* 

F i 

S ij 

S 33 

V 


A 3 

m 


F i 

/ik 

S jk 

S kk 

m 



to* 

p k _ 


where 




Vi 


4Am 


c c 

±JL 


am z A 

12 


(3.4.31) 


, . b i \ * 

ik 4 Am 


am 2 Am 

12 ~ 


V mV 

j k 


am 2 Ho Am 


The addition of all of the elements leads to a set of simultaneous 
equations written in matrix, form 

CsU [A] - [f] (3.4.32) 

where 

8 » N x N global EVP coefficient matrix 

A • N x 1 nodal EVP vector 
F ■ N x 1 nodal excitation vector 
N ■ total No. of nodes in R and on C the boundary 
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V i 

K ' 


The development of the EVP method impllcity specifies homogeneous 
Newman boundary conditions on the contour C. It remains to explicity 
specify any known nodal values of A, This is performed for an arbitrary 
node i by entering zeros into the i th row of S, except for one in the 
diagonal position and place the known nodal value of A in rhc i"“ row 
of F. 

After incorporating the boundary conditions, Equation (4.32) may be 
solved for by Gauss elimination or by an iteration scheme. The Gauss 
elimination method may be accomplished by an LU decomposition which 
also utilizes the symmetry of the banded S matrix. 

After the EVP a vector is determined, then these nodal values may 
be used in the computation of the following quantities: 

(a) H^, the elemental magnetic field Intensity which like A 

varies linearly over each element. The magnetic field intensity 
in each element may be approximated by 


H 


H + 
mo 


A i * A i + \ 


a 


(b) The element current density components in element m 
can be expressed as 

J „ ■ (c i A i + V S V' 24 ” 

J y ’ ‘ (b i A i + b ] A j + b k V /24m 

where Am ■ elemental area 

(c) The elemental eddy-current loss is 


APme 


/V ♦ j A 


and 


the total eddy-current loss is given by 
M 

Pe - l APme 
m*l 
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IV. EDDY-CURRENT DAMPING OPTIMIZATION 


INTRODUCTION 

This presentation follows very closely the work of Mikulinr.ky and 
Shtxiku^n (1) and fills in some of the mathematical gaps between their 
equations. It essentially studies the optimization of l ddy-current 
damping device containing a metal disk moving in a magnetic field of 
cylindrical symmetry. Analytical equations for the damping which is 
produced by perranent magnets for a wide range of geometrical parameters 
of this device are presented. The geometry which produces the maximum 
damping under size constraints is also obtained. 

Eddy-current dampers have numerous applications, e.g., in balances 
and in electrical supply meters. Eddy-current damping is obtained by the 
motion of a metallic body (which is tc be damped) with velocity, in a 
magnetic field. This motion creates electric current. Additionally, 
heat (Q) is produced which causes the decay of the electric current 
den&lty, J. The amount of heat produced Is given by Joule 9 8 Law, 

Q » ^ / J 2 dV* « fV 2 Watts , (4.1) 

where o is the conductivity of the body, V is the volume, an' 4 f is the 
viscous damping coefficient which we want to calculate. This coefficient 
was calculated for simple geometrical configurations by Davis and Reitz (2) 
and Schieber (3-5). A more realistic configuration, close to that used 
in some high speed levitated flywheels (6), will be studied in this paper. 

The damper is constructed of a copper disk moving in a magnetic field. 
Each of the two permanent, Identical cylindrical magnets, AA’ and BB' 

(see Figure 1), consists of two magnetic cylindrical ring3 A and A* (B and 
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B'), with radii R 2 and R^, magnetized in opposite directions. The copper 
disk, with internal radius R^ and external radius R^ can move in a 
direction perpendicular to the axis of cylindrical symmetry, C. This 
is tne motion which will be damped. 

For simplification of calculations, the hole of the disk will be 
neglected (Ro * 0) . This results in a small error which will be cal- 
culated by comparisons in two particular cases. 

Case 1 

The damping coefficient, f^» for R^ * R^ will be compared with the 
damping coefficient, f^, for R^ ■*> “ (R 3 /R 4 ■+• 0)» Both of the coefficients, 
f, and f,,, are calculated for no holes in the disk or magnets (Rq * R^ ■ 0) 
and for thin systems « R^) . For the value of R^/R^ which maximizes f, 
the difference between f ^ and f a is approximately 10Z-12%. Therefore, the 
portion of the disk outside of the magnetic field contributes no more than 
12Z of f. The part of the disk outside of the magnetic field where the 
radius of disk hole is not larger than the magnet’s inner radius (Rq <. R^) 
should contribute even less than 10Z-12Z to the damping coefficient, f. 

Case 2 

The system is thin (1^ « R^) and the magnet hole radius is large 
(R^ » R^'Rj) this case. Also, the two magnetic cylindrical rings, A 
and A* (B and B'), are equivalent in width (R^ * (R^ + R^)/2). There is 
no difference between the value of f w , calculated with Rq* 0 and R^ >*, 

and the value of , calculated for R < R^ and for R^ > R.. The hole 
in the disk still has no significant effect on the value of f. In further 
considerations, the hole in the disk will be neglected (Rq“ 0 ). 



This analysis is based on thin systems only. However, the geometry 

of real systems approaches a thin system and these approximations should 

have a wide range of validity in practice. 

A four-step procedure will be utilized to derive the damping coeffi- 

cient. First, the magnetic field, H, created by the permanent magnets A 

and B will be calculated. Using this field* the eddy-currents for small 
-*■ 

velocities, V, will be calculated, neglecting the skin effect. Calcula- 
tion of the heat production and the damping coefficient, according to 
Equation 1, is the third step. Finally, the geometry of the damper will 
be optimized *“o obtain maximum damping. 
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Reviewing some of the terns that will be used: 




Symbol 

Q 

-► 

V 

V 
f 

e 

s- 

H 

-»■ 

B 

H 

-*■ 

E 

+ 

J 

-*• 

D 

P 

* 

y 

y o 



Heat Production 

Velocity 

Volume 


Viscous Damping Coefficient 
Conductivity 

Magnetic Field Intensity 

Magnetic Induction or 
Magnetic Flux Density 

Magnetization 

Magnetic Potential 

Electric Field Intensity 

Electric Current Density 

Electric Flux Density 

Electric Charge Density 

Electric Potential 

Permeability 


Free Space Permeability 
4ir x 10" 7 


Unit 

Vatts 

Meter/ Second 
Meter 3 

Newton-Second /Meter 

(Ohm-Meter)' 1 

Ampere/Meter 

_ , Newton-Second 

Tesla ■ 

Coulomb-Meter 

Tesla 

Ampere/Meter 
Volt/Meter 
Ampere/Meter 2 
Coulomb/Meter 2 
Coulomb /Meter ^ 
Volt/Meter 
Henry's /Meter 
Henry' s/Meter 


Note : 1 Volt ■ l Joule/Coulomb 

1 Ampere - 1 Coulomb/ Second 

1 Joule - 1 Newton-Meter 

1 Ohm - 1 Volt/Ampere 

1 Henry ■ 1 Volt-Second/Ampere 
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4.1 Magnetic Field 


A magnetic field caused by stacic magnets is calculated by using 2 
of Maxwell 1 s equations* 

4 4- 4 4 4 4 4 4 n 

A*B * 0, or, since B ■ pH and H ■ -ViJj, V*Vif; * V z iJ; « - 0 

4 

-4 -4 -4 3D 

V x H * J + ~ “ 0 

a t 


(4.1.1) 


To facilitate calculations of the magnetic field, the cylindrical 
coord inate system will be used. The central axis C will be the z axis, 
r the distance from the z axis, and 0 the angle measured from an arbi- 
trary fixed direction in the plane perpendicular to the axis. 

The boundary conditions are as follows: 

On the iron yoke, s « + ij/(r,0, + ij) ■ 0 

Assume also that at r * R^, R^, 0, z) « 0 

The last boundary condition allows use of the same orthogonal system of 
functions for both magnetic potential, and electric potential, 
arriving at a simple analytic solution for f. This last boundary condi- 
tion will result in negligible error for a very thin system, « R 3 . 

The error will be small for a disk with large radius, R 4 » R 3 . It will 
be assumed that in general this error will be negligible. 

The other boundary conditions ensure that to and B are continuous on 

z 

the surface at z * + ^ 2 * 

Due to the cylindrical symmetry of the system, the magnetic potential, 
\ p, is not dependent on the angle 6 . Using the boundary conditions listed 
above, the solution of Equation (4.1.1) is the following: 


* 



C(n) Binh Bn(&l-ft2) 

Vl 2(& n R 4> 8inh *n*l 


J 0 ( ®n r)sinh &n z 


(4.1.2) 
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(4.1.3) 

(4.1.4) 


Solving for Equation (4.1.4), * 2.405, 5.520, and 8.654 would be 

the first three solutions of this equation. Therefore, 6^ * 2.405/R^, 

f $2 “ 5.520/R^, and 6 3 ■ 8.654/R 4 would be the first three values of 

2 . These values can be used in Equations (4.1.2) and (4.1.3) to find the 
potential, which defines the magnetic field inside the copper disk. 






4.2 Eddy-Current 


The electric current density, J, created by the motion of the disk 
in the magnetic field, H, is given by 

J - a [I + W Q (V x H)] (4.2.1) 

with 

- 4 - 4-4 4 

curl E * V x E ■ 3B/3t * 0, from Maxwell's equation when 

B is constant over time (therefore, E * -7$), and div J ■ <7*J ■ 0, (4.2.2) 

from the continuity equation when p is constant over time 


Taking the divergence of Equation (4.2.1), 


-+■ 

7*J 


Aq> 


7*o[E + y Q (V x H)] - 0 


- o [>E + y D 0(V x H)J] - 0 

- a [V*(-V4>) + y [$• (V xH) J] - 0 


7 2 * “ y Q [V* (V x H)H 

--4 -4 -4 -4 <4 -4 _ 

- y Q [H*(7 x V) - V«(7 x H)] 


- -y 0 V* (7 x H) 

for small velocities that are approximately constant. 

Using Maxwell’s equation with electric fields varying slowly over 

r 4 4 4 4 4 4 -4 4 

time, V y H - J and multiplying both sides by -y 0 V, -y Q V* (V x H) • y Q V*J. 


(4.2.3) 
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Equations (4.2.1) and (4.2.3) are substituted into this equation 
obtaining 

A* - -m 0 V-<oCe + H 0 (V x H)]) 

- Po<*V’C-(-V4»)-iio(V x H)D 
“ UodlV-V* - V- y Q (V x H)] 

- yoaCv.V* - uoH* (V x V)] 

” p G a V-V* 

For very small velocities, A$ ■ 0. 


Velocity is estimated as Va£«i>, where £ is the characteristic 
size of the disk and url/t is the inverse characteristic time, x» 
Equation (2.5) can be used if £«(yo«) "%..££» skin length, or 
V«(n 0 a£) _1 . When a - 10® (Ohm-M)" 1 and £-(10- 2 - lO-^m, a velo- 
city, V, << (0.1 - l)m/sec, allows use of Equation (4.2.5). 

At the disk surface, the boundary condition is 


where J is normal to the surface component of current density, j. 


The solution of Equation 4.2.5, using the boundary condition given 
in Equation 4.2.6 is 


-4 MV 


CD 

£ 


slnji 


0 


In 


R 4 h*0 sinh0. ^ (R, 2 0? -1) * , l ,p ln"4 

in l 4 In 


v) 

JaTT) COSh P in Z 


c(m) 


sinh 6 (£,-£») sinh 0 
n i z 


A£ 
n 2 


ft t 2 /fl „ v sinh 0 £. 

m*0 J. (0 R) m 1 

1 tn 4 


2 2 
"in A 


V<W 


The set 0, is given bv 
in 

V ( W * 0 


Where the Bessel function 

V^inV - * “ VemVl - 0 


(4.2.4) 

(4.2.5) 


(4.2.6) 


(4.2.7) 


(4.2.8) 
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4.3 Damping Coefficients 
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The damping coefficient, f, is found using Equations 4.1, 4.1.2, 4.2.1 
and 4.2.7. 

From Equation (4.1), 

Q . f V 2 . i / j2dV 


but using Equation (4.2.1), 


J - oCe + p (V x H)D 


(4.1) 


(4.2.1) 


therefore, 


q « fV 2 - i / cr 2 [E + p D (V x H)3 2 dV 


Equation (4.2.2; gives 


! 


: i 

. i 


t 


■+ + + 
V x E - 0, E - -V* 


This results in 


Q « fV 2 - 0 /C-V* + p Q (V x H) J 2 dV 

- o/C(V$) 2 - 2p 0 V<KV x H) + y 0 2 (V x H) 2 ]d*v 


Solving for f, 


f - i.M 2 0R lt 3 


te(n) 8±nh 

n sinh 


-161 

n 


IV* 


X ln coth v x i n AI " /2 ) 
*ln 2 -l 


hill 


3 sinh x AL + x AL 
;1 n 


sinh x L si: 
2c (m) m 


m sinh x L, 

m x 


where x^ and x^ n are the roots of the following 2 equations: 


(4.2.2) 


(4.3.1) 


(4.3.2) 


inh(x 1/2) J/(x )“1*1 
m i m I 

X. 2 x 2 jTiTxT I 
In - m 1 m J( 


- k • *. . 

* . * * 

; •* *• i , 

- i • 

k •• ‘ , ,♦ 

X 

I I *. 1 ■ 




ft? 


VV ’ °’ V (x ln> ■ W - W - 0 


(4.3.3) 
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c(n) 


2X 2 J 1 ^ X n X 2^ “ X 1 J 1 ^ X nV “ X 3 J 1 (x nV 


l - (^-v^ 


l, /R 
1 4 


(4.3.4) 


AL - Ai/R 


R./R 
i 4 


Three particular cases will now be Investigated. 

(1) For a very thin system, ■ Jl 1 /R^«l, Equation (4.3.2) results in 

*" 1 T x J X (x ) T 2 


f - 4 ttM 2 oR. 3 (—) AL 
4 L x 


- c(n) 2 

' 2 . 2, « - 21 2 . 

n x n J, (x fl ) n x ln - 1 


X J, (x ) 
I c(m) -\- 2 - 

x ln~ X m J 2< X m> 


(4.3.5) 


(2) If the radius of the disk is approaching infinity, R^-*-*, the 
parameters 

A6 X n + l~ X n and A6 ■ X ln + 1 X ln 
n R4 ln IU» 

are small, allowing replacement of the sums in Equation (4.3.5) with integrals. 
Assume that the major contribution to f arises from x , x. »1, which will 
be verfied below. 

For the Bessel functions, use the asymptotic form 


V 2) ■ “■ ‘'- 1 ! -!> 


(4.3.6) 


resulting in 


In 


£ir + im 


(4.3.7) 


J, ' (x ) a sinirm « 0 

1 m 
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The only nonzero term in the second half of Equation 4.3.5 is 

the term with x, * x because 
In m 


lim 

x m**l n x, o — x 


V 


x ln 2 ‘V 


TT- < - 1 >° 

1 x Ln 1 " 2»* ln T 


is finite. Substituting this result into Equation 4.3.5 and replacing the 
sum with the integral. Equation 4.3.5 becomes 


irM 2 q 

2 


r M ccs) 


t 


sinhB(l^ - 
sinbftf „ 


V 


(sinhBAX. + BAJt) 


(4.3.8) 


where C(0) is given in Equation 4.1.3. Therefore, as 6 goes to zero, the 
integrand in Equation 4.3.8 also goes to zero and the primary contribution 
to the integral is from the finite 8. This results in our earlier 
assumption, x n “8®4»l* 

(3) The final case to be considered is for a very thin system, 
*l/ R 4«l, v * t * 1 a disk infinite radius, R^-*°. Substituting into 

Equation 4.3.8, 


2 < 4 i ~ 

f - ttM 2 o — - — _ 

*i 2 


i 2 ) 2 


AiJ^ dB C(B) 2 


(4.3.9) 


Using the normalization equation, 

r J (8r) J (8r ') BdS - 7 ff(r - r'). 
o xn m r 

and Equation 4.1.3, the integral in Equation 4.3.9 can be calculated, 
arriving at the final analytic result 


f 


2 


<*1 - * 2) 2 



Ai (R3 2 - R l 2 ) 


(4.3.10) 
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4.4 Optimization of Damping Coefficient 


The results of the numerical calculations and optimization of the 
damping coefficient, f, will be presented in this section. The case 
that will be studied is one in which the magnets lie in a thin system, 
L^«l, contain nr hole, Rq * = 0, and are uniformly magnetized in 

one direction, Rj ” 0. The reduced damping coefficient 


♦ ■ h2 ° f-n>— ”i Ll [m2] 


(4.4.1) 


is calculated as a function of R3/R4, magnet radius/metal dial. radius, 
with R3 - 8 x 10” , according to Equation 4.3.5. The summation indexes n 
and m were varied from one to 10, giving a 10 x 10 matrix, resulting in 
a maximum at R3/R4 “ 0.3. For R3/R4 < 0.3, the calculations are incorrect 
due to large values of n and m contributing to f. when is large. For 
R3/R4 > 03, the. accuracy of this calculation is reasonably good with 
the 10 x 10 matrix. As R3/R4 decreases In this region, $ Increases. 

Tor R3 > R.1. x ■ 0, since the magnetic field ir. the disk is uniform and 
the magnetic flux in the disk does net change during its motion for a 
thin system. Therefore, it would be expected that £ -*■ 0 for R3/R4 * 1. 

It is found, from numerical calculations, that the function $ is sensi- 
tive to the order of matrix m x n at S3/S4 ■ 1. The order of matrix was 
varied frem 1 x 1 to 12 x 12. The magnitude of $ mouocomically decreases 
6 times in this interval when the order of the matrix increases. Extra- 
polating f c- s large order matrix results in $ -*• 0, which agrees with the 
prediction made above. 

As a numerical example, the damping coefficient, f ; will 1>« calculated 
for a copper disk moving between two barium ferrite magnets with 
o » 0.6 x 10“ 8 (Ohm-H)" 1 , M ■ 0.35 tesla, (A^ - A2>/A^ * 0.66, 
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LI ■ 0.8 x 10“ 2 m, and Rj » 8 x 10" 2 m. For R 3 /R 4 * 0.3, f ■ 250 Newton- 
sec/meter. If R 4 -► *, f ■ 250 Newton-sec/meter, according to Equation 4.3.10. 

Calculations of $ as a function of (R 3 /R 4 ) 2 , with R^ - 0, L^«l, 

R 3 ■ 8 x 10“ 2 m, and three values of R 3 /R 4 ( 1 , 0.8, 0.86), demonstrate that 
the maximum value of the damping coefficient, f * 220 Newton-sec/me*" <»r, 
is obtained when (R 2 /R 3) 2 * 0.5. This implies that for optimal results, 
the area of the oppositely magnetised rings should be the same. The 
magnitude oi $ , 220 Newcon- sec /meter, for the case in which Rj * R^ is 
only 10-12* less thar. the value of $ for the case of R 3 /R 4 -*■ 0, wit*, all 
other parameter;! unchanged. Therefore, increasing the disk/magnet radius 
ratio (Po/R^) is of little value. It can also be verified that th® 
assumption made earlier, that the results obtained with the disk hole 
radius (Rq) magnet radius (R^) are not considerably different from those 
results where Rq * 0. Actually, if the disk volume contained between R 3 
and R 4 changes the value cf f 3e»s than 10 - 12 %, for io -*■ <*>. it Is natural 
that the volume of the disk with radius less than R^ would change the 
value of f even less than 10 %. 

The optimal width of the copper disk, A^, can be fomd from Equation 
3.5 for thin systems. For fixed total width, of the system, the given gap 
between magnet and disk, &3 * *2 “ Ai/2, and for real systems “ 0-1 cm and 
Aij, ■ (£^ - £2)n. “ j (^1 ” ^ 3 ) • The maximum value of the damping 
coefficient, f , with respect to A£ is 

ID 


f a 

m 27 


32 J2. 

TTM OR. 


^ 1 ~V_ £ c(n ) 2 _ 

g 2 n J^fx ) 

*1 L x n J i u n ; 


n 2 i 
x* -1 
In 


m C<,t) 2 


m 


j'<x j 

ID 


- 2 . 2 / * 

x ln~ x n J l (x m>. 


(4.4.2) 



Comparing the calculations of $ using the exact Equation 4.3.2 with 



/ 


the approximate Equation 4.3.5, it is observed that the error is less 
than 10% for values oi' 0.1. When > 0.1, the error increases. 

Therefore, when ^/R^ < 0.1, the approximate Equation 4.3.5 can be used. 

As the value of increases, the magnitude of ♦ decreases. The 

reason for this decrease is that the z component of the magnetic field 

is the primary contributor to the damping coefficient, f. This compon- 
ent of the magnetic field causes the current in the plane perpendicular 

to the z direction, which is the most important factor for damping. As the 

value of A 1 /R 4 increases, the z component of the magnetic field decreases. 
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4 . 5 Conclusions 





In summary, the most general equation for the damping coefficient, 
f, is Equation 4.3.2. Calculation of Equations 4.3.3 and 4.3.3 is necessary 
for solving Equation 4.3.2. In the case where = £^/?^<<l, the simpli- 
fied Equation 4.3.5 may be used. To obtain the maximum damping coefficient, 
f, the following parameters should be used: 

- 0.5 ; A* - - l 2 « j ( 1 ^ - Si 3 ) (4.5.1) 

The first of these two parameters allows for the areas of the oppositely 
magnetized rings to be equal. The second parameter relates the thickness 
of the disk, A£, to the magnet thickness, £^ - £j>, and to the total 
thickness of the system, £^, minus the thickness of the gap between the 
magnet and the disk, £ 3 . Using the optimal parameters. Equation 4.5.1, 
the damping coefficient, f. can be calculated from Equation 4.3.10, 
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V. PROPERTIES RELATING TO EDDY-CURRENT DAMPERS OPERATING IN CRYOGENIC 
ENVIRONMENTS 

5.1 Resistivity Values of Damper Materials 

Prom the simplified basic equation for the damping coefficient 

r - B 2 1A 
d p » 

it is obvious that is maximized by a minimum in the resistivity , p . 
Thus, it is imperative that the damper material have the lowest possi- 
ble resistivity, along with the required mechanical properties for 
operation at 27 °K. 

The electrical resistivity of most pure metallic elements at 
ordinary and moderately lov temperatures is approximately propor- 
tional to the absolute temperature. It is postulated that the micro- 
scopic mechanism responsible for the temperature dependence is the 
interference to the flow of electrons caused by the thermal agitation 
of the crystal lattice. At very low temperatures, however, the resis- 
tivity approaches a residual value almost independent of temperature. 
This residual resistance is attributed to lattice imperfections and 
impurities. A small impurity has the effect of adding a temperature- 
independent increment to the resistivity. 

Alloys, as a rule, have resistivities much higher than those of 
their constituent elements and resistance-temperature coefficients 
that are quite low. For example, the alloy, 60 parts copper, 40 nickel 
(constantan) , has a room-temperature resistivity of about 44 micro- 
ohm cm while copper and nickel separately have resistivities of 1.7 
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and 6 micro-ohm respectively. Also, while the residual resistances 
of the pure metallic elements at very low temperatures are very small, 
that of constantan is about 95 percent of the room-temperature value. 

Table 5.1 shows the temperature dependence of resistance of 
several possible candidate elements for use as an eddy-current damper. 


TABLE 5.1 

Effect of Temperature on the Electrical Resistance 
of Several Pure Elements Cl] 


(Values are given as R/Rq, where R is the resistance of a 
specimen at the indicated temperature and R Q is its 
resistance at 0°C or 273°K) 


V \. Material 

Temp. Nv K 

°K 

R/Ro 

Al 

Cu 

Mg 

Ni 

Pb 

Zn 

193 

-80 

0.641 

0.649 

0.674 

0.605 

0.683 

0.678 

173 

-100 

.552 

.557 

.590 

.518 

.606 

.597 

153 

-120 

.464 

.465 

.505 

.437 

.530 

.516 

133 

-140 

.377 

.373 

.419 

.361 

.455 

.435 

113 

-160 

.289 

.286 

.332 

.287 

.380 

.353 

93 

-180 

.202 

.201 

.244 

.217 

.306 

.271 

73 

-200 

.120 

.117 


.156 

.232 

.188 

53 

-220 

.071 

.047 


.112 

.157 

.108 

33 

-240 

.049 

.012 


.089 

.075 

.041 

20 

-253 

.0427 

.00629 


.085 

.0303 

.014 


Silver has a resistivity-ratio at 27 °K of only about twice that of 
copper, but pare silver has very poor mechanical and machining proper- 
ties. Copper has the lowest resistivity at the required working temper- 
ature and itb mechanical characteristics are satisfactory, although it's 
machineability is poor. 
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A literature search for values of resistivity, p, for a "pure" 
copper at low temperatures has indicated the values as shown in 
Table 5.2. 


TABLE 5.2 

Resistivity of "Pure" Copper 


Source 

T 

rm 

20 6 C - 293°K 

67 °F - 527 °R 

T LN 2 

-197 °C * 76°K 

-323 °F - 137°R 

Tssme 

-247°C - 27°K 
-411°F - 49°R 

A- 

...» 

_9 

-2 x 10 n-m 

-3 x 10 fl-m 

B 

— 

— 

-3 x 10- 11 n-m 

C 

1.553 x 10- 8 n-m 

1.86 x 10“9 n-a 

-4 x 10 -11 n-m 


(+.005) (273®K) 

(+.02) 


D 

1.7 x 10“ 8 n-m 

*3 x 10~9 

< 1 x 10 - ^ n-m 

E 

1.55 x 10" 8 fl-m 

-2 x 10“9 n-m 

-1.6 x 10 - * 8 n-m 

F 

-1.6 x 10“ 8 n-m 

~2 x 10“^ n-m 

-1.4 x 10' 10 n-m 


(A) L. A. Hall. "Survey of Electrical Resistivity Measurements on 
16 Pure Metals, in the Temperature Range 0 to 273°K", N.S.B. 
Technical Note 365, Washington, D.C. 

(B) P. K. Moussouros and J. F. Kos. "Temperature Dependence of the 
Electrical Resistivity of Copper at Low Temperatures", Can. J. 
Phys., Vol. 55, No. 23, 1977, pp 2071-2079. 

(C) F. R. Fickett. "A Preliminary Investigation of the Behavior of 
High Purity Copper in High Magnetic Fields", N.B.S., Cryogenics 
Division, U-235, June 1972. 

(D) "Handbook of Thermophysical Properties of Solid Materials", 

Vol. I, 1961. 

(E) D. L. Grigsby. "Electrical Properties of Copper, Manganln, 
Evanohm, Cupron and Constantan at Cryogenic Temperatures", 
Hughes Aircraft Co., October 1966, Electronic Properties 
Information Center, I.R. No. 40. 

(F) R. Barron. "Cryogenic Systeas", McGraw-Hill, 1966. 


5-3 



5 . 2 Behavior of High Purity Copper in High Magnetic Fields^ 

At low temperatures, the resistivity of copper increases 
almost linearly with increasing magnetic field. This if true 
for a range of purity of 200 < RRR < 7000 and for a tempera- 
ture range cf 4K < T < 35K. RRR stands for the Residual Resis- 
tance Ratio which is equal to R(273K)/R(4K) and is a sensitive 
indicator of purity, i.e., increasing ratio represents an in- 
creasing purity. Increases in the resistivity by a factor of 
120 have been observed for a very high purity copper sample at 
4K in a field of -100K gauss or 10 tesla. 

However, for any magnetic field and temperature, the value 
of the resistivity can be accurately predicted by Kohler’s 
rule which is 



where B is the magnetic field (flux density), R q is the resistance 
at zero field and AR * R(B)-R o « For pure metals, f is a single- 
valued and monotonically increasing function of B/R q . 

A Kohler diagram for copper is shown in Figure 1, which is a 
plot of AR/R^ versus B*RR(T) , where RR(T) is the resistance ratio, 

R(273K)/R(T), and R(T) - R . 

o 


1. F. R. Fickett. "A Preliminary Investigation of the Behavior of 
High Purity Copper in High Magnetic Fields", Cryogenics 

Division, U-235, June 1972. 
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Consider now an example for a practical damper design with 

a copper disk. Using a realistic magnetic flux density for a 

4 

permanent magnet of one tesla or 10 gauss, an extremely low 
value of resistivity at 27K = 4 x 10 ^ 0~m and a resistivity 
at 273K - 1.55 x 10~ 8 0-m, then B*R(273K)/R(T) « 3.9 x 10 3 
kilogauss and from Figure 1, we see that the resistance due to 
the magnetic field has increased by 100 percent or a factor of 
2. However, if a more conservative value of resistivity at 
27K of 5 x 10”^ 8 0-m is used, then the resistance increase is 
only about 3 to A percent. 


FIGURE 1 

Kohler Diagram for Copper 


A R 

K 



B-R(273K)/R(T) - (K gauss) 
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Thus, it would appear that the dominating factor in the 
increase of resistivity due to the magnetic field would be 
the purity (RRR) of the copper • A practical value of resis- 
tivity at 27K would probably lie somewhere between the two 
examples above and thus the problem of an increase in re- 
sistivity due to a magnetic field and hence a decrease in 
the damping coefficient should not be a severe or overriding 
concern. 
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5.3 Eddy-Current Depth of Penetration 






1,2 


i t 

i 

; i 


* ' 


The solution of Maxwell's equation, when a sinusoidal magne- 
motive force is applied on a nonferrous conducting plate, gives a 
wave of flux density that enters the plate from the outside 

surface at the start of each half-cycle, and penetrates to a depth 6, 
called the depth of penetration. If 6 is less than d, the half 
thickness of the plate, the magnetic flux as well as the eddy- 
currents generated are essentially restricted to a layer of depth 
6 on each surface of the plate. If the 6 is larger than d, the 
flux-density waves from each side meet in the center of the plate 
before the end of the half-cycle of the sinusoidal magnemotive 
force and eddy-currents flow throughout the full width of the piste. 

The equations listed below, 

7 2 h 


7?E 


- jwopH 
* jdJOyE 
V 2 J - jwcrpi, 

give the basic relation between time and space deri. tives of the 
magnetic field, electric field, or current density for any point 
located in a conductor. 

Solving for the current distribution equation, for the case of a 
plate conductor vith current flow parallel to the surface, the depth 
of penetration is 

S meters. 



1. P.D. Agarwal. "Eddy-Current Losses in Solid and Laminated 
Iron", AIEE Trans., 1959, 78, Pt. (1), p t > 169-181. 

2. S. Ramo and J. Whinnery. "Fields and Waves in Modern Radio", 
John Wiley & Sons, Inc., NY, 1953. 


' . ill, 
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The complete solution indicates that the current magnitude 
decreases exponentially as it penetrates into the conductor. Thus, 

6 is the depth for which the current density has decreased to 
1/e (~36.9) of its value at the surface. 

From the standpoint of an eddy-current damper, it is apparent 
that there would be an optimum range of thickness for the conducting 
plate, i.e. , too thin a plate for the frequencies expected and the 
full damping potential is not utilized, and with too thick a plate, 
there is excess plate material and weight, etc., which is serving 
no useful purpose. Using a conservative value of p ■ 0.5 x 10 Q-m 
for copper at 27 °K and y * 4tt x 10 ^ Henrys/meter , a plot of the 
depth of penetration versus frequency is shown in Figure 1. 

This figure shows that for the range of frequencies of interest 
for the SSME, namely, ^36,000 RPM or 600 Hz, the depth of penetration 
is only on the order of 5 x 10" 2 cm. Even at 1000 RPM or 16 Hz, the 
penetration depth is quite small, 6 * 0.15 cm. Thus, at operating 
speeds of ~36,000 RPM, with aii 0.25 inch (0.635 cm) copper disk as 
the damping conductor, l >st of the eddy-currents would reside very 
the surface and little use would be made of the bulk of the material. 
It would seem that, at first glance, perhaps a much more efficient 
design would consist of a layered or ^aminated type of disk con- 
struction, such as shown in Figure 2. 


FIGURE 2 

Laminated Damper Disk for Improved Efficiency 
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Several thin copper sheets, with insulated layers between, could 
make up the total disk thickness, with the thickness of one copper 
sheet still much greater than the depth of penetration, hut now the 
total magnitude of eddy-currents has increased by a factor equal to 

i 

the number of laminations of sheets. This, in turn, would mean that 
the damping effect has also increased by the same factor. 

] 

This design is exactly opposite to the design of transformer j 

cores, etc., where the laminations are made parallel to the changing \ 

magnetic field in an effort to eliminate or drastically reduce the j 

eddy-current and hysteresis losses. j 
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5 . 4 Magnetic Induction Variation With Temperature 
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The remanance, or magnetic induction which remains in a mag- 
netic circuit after an applied magnetomotive force is removed, is 
temperature dependent. Generally, it will decrease as the tempera- 
ture increases and will become zero at the Curie point, at which all 
ferromagnetic properties vanish. This, then, implies that the 
remanance would conversely increase as the temperature decreases, 
and indeed this effect is observed in some cases. There are both 
non-reversible variations and reversible variations as a function of 
temperature. 

The non-reversible effect results in a change in the remanance 
of a magnetized magnet and it's circuit which has been temperature 
cycled. This non-reversible change on stabilizing processes is 
associated with a loss in the remanant induction. But the initial 
value of remanance may be restored by remagnetizing the stabilized 
magnet, so long as the temperature variation did not result in an 
Irreversible metallurgical change of the magnet material. 

After stabilization over a given temperature range, any further 
changes of remanance within this temperature range are reversible. 

This relative reversible variation in percent is calculated by measur- 
ing the remanance B^(t) at the temperature to within the stabilized 
range and then comparing this with the room temperature remanance 
B^(20) by the expression 


B d (t) ' B d <20) 
B d (2C) 


x 130% 


Thus, a temperature coefficient over a given range can be determined. 
Both the non-reversible and reversible variations are dependent upon 
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a shape factor or L/D ratio (length of magnet to the equivalent 
diameter of magnet) and to the magnetic material type. For 
instance, for a temperature range from -60°C to +C0°C, the net 
coefficient is negative and may be as large as several tenths of 
a Z per °C. Alirsco V, for example, has a temperature coefficient 
of -0.024 %/°C for L/D * 8 over a temperature range of 0 to 80°C. 

Mr. Shuk P.ashidi of Hitachi Magnetics Corporation indicated 
that they test magnets down to -100°C and he estimated that the tem- 
perature coefficient may be as large as 0.04 Z/°C at very low temper- 
atures. Their catalog gives a value of 0.033 Z/°C for HICOREX, a 
rare-earth cobalt permanent magnet, from 20°C to -100°C. 

Thus, while hard data on magnet properties at low temperatures 
is apparently not presently available, it appears as though a bene- 
ficial Increase of magnetic induction and hence the damping coefficient 
may be experienced. For example, with a cobalt magnet and increase 
in induction of approximately 10Z or more should occur at 27 °K. 

Since the damping coefficient varies as the square of the magnetic 
induction, the net result should be at least a 20Z increase in 
damping. Laboratory testing of this phenomena should be carried out 
to establish more accurate data on the temperature coefficient and 
the effect of temperature recycling on the magnetic properties. 
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VI. EXPERIMENTAL TESTS 


A series of tests, using a simple vibrating rod, was conducted to 
observe the phenomena of eddy-current damping and the effect of damper 
material thickness, layers and temperature. 

6.1 Vibrating Rod Set-Up 

The basic apparatus consisted of a 1/2" thick aluminum plate, which 
was securely attached to a rigid upright stand [see Fig. 6.1]. The 
rod length was adjustable. A half -flat on the end of the rod allowed 
for the attachment of the damper. The damper was always 2" wide x 3 1/2" 
high and the thickness was varied from 0.031" to 0.275". 

A rather ancient magnetron Alnico V magnet provided the magnetic 
flux. It had a 1.25" diameter pole face with a gap of 0.82" and produced 
an average field of approximately 1200 gauss (0.12 tesla). The magnet 
was mounted on a small laboratory ack and could be raised so that it 
was centered on the damper, or J t could be lowered about 6" to make tests 
without the magnetic field acting on the damper. 

The relative amplitude of vibration of the rod was measured by a 
proximity probe mounted so as to detect motion about 2" above the damper. 
The output of the proximitor was connected to one input of an HP-5420A 
Digital Signal Analyzer. The rod was struck to induce vibrations by an 
impact hammer with an attached accelerometer. The accelerometer output 
was connected to the second input and this permitted a variety of infor- 
mation to be obtained directly rrom the HP-5420A, such as frequency 
spectrum, transfer function — both rejl and imaginery, time average for 
observing the signal decay and the damping. 
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A styrofoam container, with and without the magnet centered on 
the bottom, was filled with liquid nitrogen, placed on the jack and 
raised and lowered for observing vibration data for low temperature 
(77 k) damper tests. 

i 

l It was found that, within experimental error, the observed damping 

' effect did not change for tests, without the magnet, conducted at room 

I temperature or at LN^ temperature. 



t 



6-3 







6.2 Damping Calculations 

The damping coefficients for the vibrating rod were 
calculated from the rod and damper mass, the resonant fre- 
quency and the log decrement of the rod motion wave-form 
after being struck. The critical damping of a one- 
dimensional sy 'tern is related to the stiffness, mass, and 
natural frequency by 

C = 2 /{cm ~ 2m<i) . 
c n 

The actual or measured damping coefficient, C^, is 
related to the critical damping by the damping ratio or 
factor, 


and for small values of this damping ratio, 


where 5 


log decrement « In *1. - In x 0 , 

v n — 


x *-* amplitude of vibration, 
n * number of cycles. 

Thus, the measured damping coefficient can be determined 


from 


C m - cC - 
m c a’ 


2m(2irf n ) 


C « 2m 6 f r 


lb-ccc 


N-sec 







The theoretical damping coefficient, as derived earlier 
in this report, may be calculated from 

( 6 . 5 ) 

where B * magnetic flux density, 

V * volume of damper material in magnetic field, 
p * resistivity of damper material. 

The theory thus indicates that the damping effect should 
vary linearly with damper thickness, at least for thickness 
less than the penetration depth. Also, the damping should 
increase by a factor of about 8 when the damper material is 
at liquid nitrogen temperature, since pT^ c 1.6 x 10” 8 0-m 
and pT TN * 2 x 10 -9 n-m. 






6 . 3 Practical Observations 


The exact nature and paths of the eddy-currents are 
extremely difficult to examine analytically, as several of 
the texts on motors and E & M have acknowledged. 

For the case of a solenoid or iron- core transformer , 
the eddy-current picture is fairly clear. The changing 
magnetic field and the volume of material in which the eddy- 
currents are induced are both well definc-i and thus the 
theoretical damping coefficient relation, Eq. 6.5 ; should 
yield accurate information. 

However, for the proposed SSME eddy-current damper con- 
figuration, the picture is considerably different and unclear. 

The EMF which produces the eddy-currents is induced solely by 
a changing magnetic field, according to Faraday's Law: 

EMF - MX . (6.6) 

3 1 3 t 

Thus, for the case, where the damper material extends in 
all directions well beyond the area of the magnet pole faces, 
it is obvious that the magnetic field is not going to be 
changing in the moving material over most of the internal 
area of the pole faces, since the field is assumed uniform 
here. The magnetic field will be changing in the moving 
damper material near the edges and fringe regions outside the 
pole face area where the stationary magnetic field exhibits a 
gradient or non-uniformity. More simply expressed, if there 
was no fringing and the magnetic field was uniform exactly 
over the entire pole face area, then EMFs (hence eddy-currents) 
would be produced in the moving material only at the edges of 
the pole face area where change of flux is occurring, and no 



EMFs would be generated in the material interior area of 

the pole face, since there is no in this region. The 

implications of this indicate that it would be extremely 
difficult to calculate a theoretical damping coefficient from 
Eq. 6.5 for this practical situation because, due to fringing, 
neither the magnetic field B nor the exact damper material 
volume (area x thickness) is known for the regions of inter- 
est. 

To check the validity of the concept that no eddy- 
currents would be produced in the moving damper material in 
the interior area of the magnet pole face, a simple qualita- 
tive type experiment was performed. A one foot long insulated 
pendulum, with a circular 1/32 inch thick copper disc of dia- 
meter 5/8 inch attached at the bottom, was rigged so that it 
would swing for 65 to 70 cycles in air when started with an 
initial deflection of 1/4 inch (see Fig. 6.2). 

Next, a magnet with a 1 7/16 inch diameter pole face was 
placed so that the copper disc on the pendulum was halfway 
into the pole face area and would swing in the center of the 
gap separation perpendicular to the magnetic field. Thus, a 
portion of the disc would be cutting magnetic lines as it 
travelled from a low to high magnetic field and vice-versa, 
and an eddy-current damping effect should be present. Indeed, 
this was observed as the pendulum would now swing for only 
35 to 40 cycles with an initial deflection of 1/4 inch from a 
rest position. 
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Finally, the magnet was again moved such that the disc 
was now in the center of the magnetic field. The swing test 
was then found to again produce 65 to 70 cycles for the 
1/4 inch initial deflection. These 3 cases were each repeated 
for about 10 trials. The results are shown in the table below. 


TABLE 6.1 

Pendulum Swing Test Results 


Initial Conditions 

No. of Cycles of Swing 
for 1/4" Initial Deflection 

No Magnet - Free Swing 

65 to 70 

Copper disc on pendulum half 
into magnet gap - large 
gradient present 

35 to 40 

Copper disc on pendulum lo- 
cated at center of magnetic 
field (i.e., gap) - uniform 
magnetic field 

65 to 70 


This brief experiment clearly demonstrated that the area 
or volume of damper material in the interior uniform magnetic 
field region does not contribute to eddy-current damping, since 
there is no variation of magnetic flux within the damper material 
to induce the EMF ' s needed for eddy-current generation . Thus , to 
optimize the damping effects for this configuration, which is 
similar to that proposed for the SSME, the regions of largest 
magnetic field gradients should be maximized, along with the 
volume of damper material in this region. The magnetic shape 
could have a strong influence on the fringing or gradient 
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picture, i.e., a square cornered magnet should produce more 
fringing as compared to a round magnet. Also, it would appear 
that several long rectangular magnets would be more effective 
as compared to a circular magnet of equal pole-face area, 
since the area of uniform magnetic field is greatly reduced, 
whereas the magnet edge area, at which the large gradients 
occur, is greatly increased. 

An idea for possible exploration would be to insert small 
segments or pieces of a ferrous material into the copper damper 
in the inner pole-face region where the flux density, B, is 
generally very uniform, and hence produces no induced emf's 
for eddy-current damping when the damper is moving. With a 
ferrous insert, the magnetic field would now have a gradient, 
thus inducing emfs which produce the desired damping effect. 

The inserts should not penetrate the total dampen thickness, 
so that some low resistivity copper is present to maximize 
the eddy-current damping. 

Again, this is an extremely difficult problem, to calculate 
the damping from theory. Probably the mos^ feasible path would 
be for a series of experimental tests to be conducted with 
various empherical relations resulting. 
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6.4 Results of Variation of Volume and Temperature of Damper 

Material 

Using the vibrating rod set-up as described in Section 6.1, 
a number of tests were conducted to observe the effect of 
eddy-current damping with variation of thickness and temperature 
of the damper material. The measured damping coefficient, 

(Eqn. 6.4), was determined from the calculated mass and the log 
decrement and vibrating frequency as observed on the time-record 
plot from the HP-542 OA. Typical plots of the time-record taken 
without and with the magnet are shown in Figs. 6.3 and 6.4. 

The increase in damping due solely to the effects of 
eddy-currents could thus be determined by noting the difference 
in the measured damping coefficients, C^, from the "without" 
magnet to the "with" magnet case. This increase in damping, 
expressed as a percent, is plotted in Fig. 6.5 as a function of 
damper thickness for both the room temperature and liquid 
nitrogen tests. The increase in damping, AC m in lb-sec/in, 
versus darnoer thickness is shown in Fig. 6.6. 

B 2 LA 

If the theoretical damping coefficient, C T = — - — , is 
compared to the measured increase in damping, AC^ = (magnet) 

- C (no magnet) , by a factor K, such that 



then the plot of K versus damper thickness as shown in Fig. 6.7 
gives an indication of how accurately this relationship 
represents the response of the set-up configuration, compared 
to an ideal one, where K should equal unity. 


6-li 

















ORIGINAL PAGE IS 
OF POOR QUALITY 


% INCREASE IN DAMPING ( with and w/o magnet ) 

vs 

DAMPER THICKNESS 


B av g = 1200 Gauss 
Copper 



insulated spacers 





(Ib-sec/ln) 


I 


ORIGINAL PAGE 19 
OF POOR QUALITY 





Ki 

i 

■f 


nr 



figure 6.6 



FACTOR K - (C t /AC m ) 




o 100 200 300 

DAMPER THICKNESS - (MILS) 

Figure 6.7 
6 - 1 6 


1 


t 

l 


1 , 


\ 


( 



i 1 * * 


r 




From an inspection of these three figures, several 
immediate observations can be made. 

1. The eddy-current damping does increase with 
increasing damper thickness, although not 
linearly as the theory predicts. 

2. There is a large increase in the eddy-current 
damping as the temperature is lowered from room 
temperature to liquid nitrogen tempera tux'e , 76 °K. 
However, for the smallest thickness, the magnitude 
of the increase was only about one-half the expected 
value of 6, and as the damper thickness increased, 
this value decreased even more. 

3. The large values of the factor K indicate that the 
measured eddy-current damping is considerably 

less than the theory predicts and that this discrep- 
ancy increases with damper thickness. 

4. The two data points. using a laminated damper, made 

of individual 0.032 inch thick copper damper sections 
with a thin paper sheet between sections, perhaps 
suggest that there may be an additional increase in 
damping achieved by taking into account the eddy- 
current depth of penetration as explained earlier. 
However, more experimental data is necessary to 
further explore this concept. 


6-17 



6.6 Discussion 


The results of this rather crude test show that the 
measured eddy-current damping does not compare favorably with 
the expected theoretical damping values for the various 
tests. There are several obvious reasons relating to the 
causes for this discrepancy and they are as follows: 

A. The proper value of the magnetic flux density, B, 
to be used in the theoretical deunping coefficient 
is very difficult to determine. For this experi- 
mental set-up, a profile of the flux density at 
the mid-span between the two pole faces showed that 
at the center, B os 1700 gauss, while around the 
mid-span circumference, B *= 1000 gauss. A value 
of 12 00 gauss was used in the calculations. However, 
as noted in Section 6.3, Practical Observations, the 
eddy currents necessary for the desired damping 
are induced by the moving damper only in the presence 
of a magnetic field gradient, i.e., in the near-regions 
to the edges of the magnet pole faces. Eddy currents 
are not induced in a constant field. Thus, with the 
fringing effects being an unknown for this set-up, 
it would be very difficult to assign a number for the 
flux density, B. With the proper equipment, a 
complete profile of the flux density at the plane 
of the damper could be obtained over an area consid- 
erably greater than just the pole face area, and 
perhaps an average or effective value of B could be 
determined. 
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B. the proper value of damper area to be used in the 
theoretical damping coefficient is also difficult 
to determine. The same reasoning as discussed in 
A above applies to the area parallel to the plane 
of the pole-face. The correct value of area should 
be that over which there is a magnetic field gradient 
occuring. Some effective value of area should be 
used, but the magnet pole face area was used for 
these tests. 

C. For the low-temperature tests conducted with liquid 
nitrogen, it is possible -hat the damper did not come 
to an equilibrium temperature, since no temperature 
measurement of the copper material was made. Within 
about five minutes after insertion into the liquid 
nitrogen, the violent boil-off from the nitrogen 
vaporization .;ad almost ceased and then data-taking 
began, but the actual temperature of the damper may 
have been higher than the temperature used in the 
calculations, due to insufficient time for thermal 
equilibrium, and thus the damper resistivity would be 
greater than value used for the calculations. 

Also, the calibration of the displacement probe is somewhat 
suspect at the low temperatures. This type of probe is 
temperature sensitive and the gap or D.C. voltage was noted 
to change from 8.2 volts at room temperature to over 11.2 volts 
at LNj temperature. The probe was approximately 1 to 1 1/2 
inches above the liouid nitrogen level and it definitely was 
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much colder than the room temperature. However, no attempt 
was made to determine its temperature or to calibrate the 
probe at low temperatures. The change in calibration for a 
70°F to 100°F change is given as about 1% by the manufacturer, 
but no data was available for the lower temperatures. 
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CRITICAL SPEED ANALYSIS OF NASA EDDY-CURRENT DAMPER TEST APPARATUS 
ackground and Introduction 

: igure 1 represents the schematic diagram of the NASA edoy- 
lt damper test apparatus. The rotor configuration is designed 
irate in liquid nitrogen. The object of the eddy-current 
• test apparatus is to examine the damping characteristics 
•assive eddy-current damper operating in cryogenic conditions. 

.\jp' t of this research program is to determine the feasibility 
>e application of such a device for cryogenic turbomachinery, 
is the liquid oxygen and hydrogen pumps used on the SSME. 
in the development of high speed high performance turbopumps 
ing with cryogenic liquids such as oxygen or hydrogen, it is 
ely difficult to incorporate damping into the bearings or 
For example, modern aircraft engines mounted on ball bearing 
ts must incorporate squeeze film oil dampers in order to control 
orations caused by unbalance response, or the nonsynchronous 
caused by aerodynamic cross coupling effects. Under cryogenic 
ons, it Is impossible to incorporate the conventional squeeze 
imper design. High performance turbopumps that run through 
it frequencies must have additional external damping Incorpor- 
to the bearings or structural system in order to adequately 
the vibrations. This is necessary if high bearing life of 
ling elements is to be assured. 
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FIGURE 1 - EDuY-CURRENT DAMPER TEST APPARATUS 


In this phase of the investigation, the critical speeds of the 
NASA eddy-current test apparatus are examined to determine the mode 
shapes aid percent of strain energy distribution between the shaft 
and the bearing supports. The eddy-current damper will be mounted 
at the No. 2 bearing location (see Fig. 2). The eddy-current damper 
will be most effective if a high percentage of the system strain energy 
is associated with the No. 2 bearing support. 



7 . 2 Critical Speed Analysis of Original Design 


The '1ASA eddy-current test apparatus was analyzed, using the 
computer program CR7SPD developed by the University of Virginia to 
operate or the HP-9845B computer system. Incorporated with the 
computer program is the graphics procedure to illustrate the mode 
shapes. Figure 2 represents the rotor model of the NASA eddy- 
current test apparatus. For the first design, the system was con- 
sidered as a two-bearing system in which the first bearing is a 
relf-aligning ball bearing with an estimated stiffness of 85,000 lb/in. 
The eddy-current damper will be located at bearing 2 and its stiffness 
will be determined by the comb, nation of the retainer spring rate 
stiffness plus the additional stiffness generated by the magnetic 
field. For this bearing, a value of 10,000 lb/in was chosen. 

The first model has 14 stations and the characteristics of fhis 
model are given in Table 1. The total weight of the test rotor is 
approximately 16.6 pounds. The system was examined for a speed 
range of 200 to 20,000 RPM. There is predicted to be only or»<? critical 
speed in the operating speed range and this value is 3,656, as noted 
in Table 2. This table represents the first mode shape and also 
gives the distribution of strain energy in the shaft md bearings 
and the distribution of rotor kinetic energy. It is of interest to 
note that in the original design, 73* of the total strain energy is 
in the strain energy of bending and only 27 % is in the bearings. 

This distribution of strain energy is not particularly desirable 
and It is perferable to have a higher percentage of strain energy 
in the damper bearing location. It is also of interest to note 
that sections 4 through 6 contain over 67% of the strain energy of 


7-4 









TABLE 1 


msfl EDDY CURRENT DP'IPFR TEST APPARATUS 
TQUIL N2 SYSTEM-1 5,600 PPM DESIGN , STIFF SHAFT 

N« 14 NBRG* 2 NCFlSE= 1 NMODES* 3 Eps* . 

I M(LBS) L ( I ) D< I ) El IP 


OR' NAL PAGE S3 
OF POOR QUALITY 

0006 1 0 
IT 


1 

. 1 

2.70 

.59 

9. 204E+04 

23. 442E-04 

46. 74E-03 

2 

. 1 


.59 

1.784E+05 

28. 997E-04 

47. 14E-03 

3 

.0 

.33 

.59 

1 . 784E+05 

11. 1 10E-04 

78. 72E-05 

4 

.5 

4.42 

1.00 

1 . 473E+06 

6 1 . 957E-03 

83. 08E-02 

5 

1.0 

4.42 

<s> 

© 

1 . 473E+66 

12. 286E-62 

16. 61E-01 

6 

.6 

.68 

1. 15 

2. 576E+06 

77. 923E-03 

84. 25E-02 

7 

6.6 

1.62 

1. 19 

2. 953E+66 

58. 303E+08 

40. 10E+00 

8 

6.7 

1.00 

1.50 

7.45SE+06 

1 1 . 546E-02 

13. 43E-02 

9 

.4 

1.49 

1.60 

1 . 473E + 06 

90. 956E-63 

96. 64E-03 

10 

.3 

1.49 

1 . 66 

1 . 4 73E «-66 

4 1 . 258E-03 

81.29E-83 

11 

.2 

.28 

.79 

5. 736E+05 

22. 144S-03 

41 * 53E-03 

12 

.0 

.28 

.79 

5. 736E+05 

30. 301E-04 

17.69E 04 

13 

.0 

.47 

. 75 

4. 659E + 05 

35.677E-04 

24. 41E-64 

14 

.0 

0.00 

.75 

4 . 659E+05 

26. 527E-04 

15.57E-04 


ROTOR WEIGHT- 16.6 LB, LENGTH* 19.5 IN., C.C.* 13.0 IN. FROM LE r T 
Ipt *5. 89E+01 LB-IN /S 2, Itcg*8. 32E+02 D*vg* .9 IN., EI*vg-l . 53E+06 
BEARING STATION LOCATIONS : 

BRG. NO. 1 - 3 BRG. NO. 2 * 12 

BRG. STIFFNESS VALUES : 


BRG. NO. 1 ST. 3 K- 85606 L.B-'IN 
BRG. NO. 2 ST. 12 - K- 10000 LB' IN 

ROTOR BOUNDARY CONDITIONS ARE! FREE-FREE 
WHIRL MOTION IS S 

Irp»« 500 Drpm- 1000 Frpm- 30000 


L1AVG- 3. 383E-06 L2AVG- 4.334E-06 L3AVG- 4.333E-06 ElftVG" 1 . 529E+06 
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TABLE 2 

ORlGiKML PAQ£ S3 

NfiSfl EDDY CURRENT DAMPER TEST fiPF'fiRATUS OF F>OOR QUALITY 

LIQUID N2 SYSTEM-13,000 RPM DESIGN RANGE 

UNDAMPED ROTOR MODE SHAPES AND tnEFG, DISTRIBUTION 
WITH TRANSVERSE SHEAR DEFORMATION 
SYNCHRONOUS FORWfiPD MODE SHAPE 


NO. 1 CRITICAL SPEED * 3656 ITER* 6 DELTA*-. 000001993 



ST 

X 

THETA 

M 

V USHfiFT UBEARING 

Kbra 

Ket 

Ke~ot 

1 





(DIM STRAIN ENERGY) 

(Dlii l 

KINETIC) 

f. 

1 

-.609 

.311 

0. 0000 

0. 0000 


0 

b 

j 

2 

-.050 

. 310 

“.0005 

-.0002 


0 

a 

) 

3 

.018 

.310 

- , 000b 

-.0002 0 1 

85,000 

0 

-? 

f 

4 

.086 

.309 

-.0058 

-.0200 


0 

d 

l 

5 

.870 

. 188 

-.0751 

-.0199 10 


3 

3 

i 

6 

1.000 

-. 133 

-. 1388 

-. 01S3 57 


3 

a 

i 

7 

.939 

-. 138 

-. 1476 

-.0169 1 


54 

-1 


8 

.783 

-. 150 

-. 1368 

.0136 3 


38 

a 


9 

.68' 

-. 152 

-.1058 

. 0396 


2 

a 

>* 

1 

10 

.5? 

-. 164 

- . 0580 

.0410 2 


1 

a 

► 

11 

. 3L8 

-. 169 

-.0092 

.0419 


0 

a 


12 

.326 

-.170 

.0000 

.0422 0 26 

10,000 

0 

-a 

r. 

13 

.294 

-. 170 

.0000 

-.0001 


0 

-a 

y 

14 

.241 

-. 170 

.0000 

-.0000 


0 

-a 

•» 

r 





7 3 2? 


101 

-l 

> 

Utot&l 

*20. 30E+02; 1 

Ke total 

= 20. 36E+82; 5s ERROR ENERGY 

BALANCE* 

-. 3 


CRITICAL SPEED SUMMARY 


NASA EDDY CURRENT DAMPER TEST APPARATUS 
LIQUID N2 SYSTEM-15,000 RPM "PSIGN RANGE 
WITH TRANSVERSE SHEAR DEFORMATION 

SYNCHRONOUS CRITICAL SPEED ANALYSIS 

Brg. NO. 1 ST. 3 K« 05000 Lb/ In 
Brg. NO. 2 ST. 12 K* 10000 Lb/In 


NO. CRITICAL SPEED Wmode Iruode UTruode Kmode Ushat't Ubr j KEt KEr 
RPM (HZ) LB LB- 1 N A 2 LE LB/'IN (DIM. STRAI 4) (DIM. KIN) 

(It-v'w*Ip> (ENERGY) 


1 3,656 ( 61 > 10.8 -.1 10.7 4,072 73 27 101 -1 





7-7 



bending. This indicates that there is considerable bending at this 
location for the first mode. The strain energy at the damper bear- 
ing may be increased by increasing the effective shaft diameter at 
locations 4 through 6. 

Examining the kinetic energy distribution, it is seen that 
only one percent of the energy is associated with gyroscopic effects. 
The rotor total kinetic energy may _ viewed mainly as translatory 
motion, rather than gyroscopic motion. This implies that a simpli- 
fied rotor modal may be used to represent the dynamic characteristics 
of tne experimental test apparatus for future analysis and calcula- 
tions. This feature is highly desirable when running the experimen- 
tal facility in order to experimentally determine the effective 
damper stiffness and damping coefficients. 

It is also of interest to "ote that 92% of the kinetic energy 
is associated with stations 7 to the large disk on the rotor. 
Therefore, this rotor may be successfully balanced for the first 
mode by only a single plane at the major disk location. In addition 
to the rotor kinetic energy, the effective rotor modal weight and 
modal stiffness is given in Table 2. 

Figure 3 represents the rotor first mode. From Table 2 and also 
Figure 3, it is noted that the maximum amplitude occurs at station 6. 
Figrre 4 is similar to Figure 3 In that it represents an animated 
mode shape for the first mode. Note that the amplitude of notion 
at the ball bearing support is almost a node point. It is also of 
interest to observe that the bearing amplitude at the damper is only 
33% of the maximum rotor arrp.'tude. This is an indication that the 
original shaft, design may be too flexible. 
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NRSR EDDY CURRENT DRMPER TEST RPPRRRTUS 
LIQUID N2 SYSTEM- 1 5 , 000 RPM DESIGN RANGE 


UNDAMPED SYNCHRONOUS 5HRFTM0DE5 
Wt- 15.4 LB Lt- 19.5 IN. 
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7.3 Stiffened Rotor Design 


A design rotor for the eddy-current damper apparatus was con- 
sidered in which the shaft section between 4 and 6 was increased 
from0.62to 1.0 inches. Figure 5 represents the stiffened rotor 
design. Table 3 gives the critical speed mode shape for the rotor 
with the stiffened design. The critical speed has been increased 
to 5,441 RPM. Note that the strain energy at the damper support 
has increased from 27 to 60%. The performance of this model would 
be much more satisfactory in the test apparatus. The stiffened 
rotor was analyzed to 30,000 RPM. Only one critical speed was 
determined to be in the operating speed range. The mode shape and 
animated mode shape are shown in Figures 6 and 7, respectively. 

A third model was run in which a coupling weighing approximately 
one pound was placed at the shaft end. A timer-pulley arrangement 
will be used to drive the rotor system. The rotor model with the 
pulley arrangement is shown in Figure 8. The analysis of the rotor 
critical speed with the pulley indicates that the overhung pulley 
will have very little effect on the rotor first critical speed. 

The balancing should be done primarily at the major disk location. 

The Influence of the pulley is to cause a second critical speed 
to occur, just outside the operating speed range at approximately 
21,000 RPM. This frequency may be changed somewhat depending upon 
the exact weight of the pulley. The shaft modes and animated modes 
are given in Figures 9, 10 and 11, respectively. 
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7.4 Summary and Conclusions 


The critical speeds of the NASA eddy-current damper test rig 
were calculated. It was determined that there would be only one 
critical speed present In the operating speed range of 15,000 RPM. 
However, in the analysis of the original rotor design, It was 
determined that the eddy-current damper bearing wc Id contain only 
27% of the strain energy of the total system for the first mode. 

This low percentage of strain energy would make the eddy-current 
damper extremely inefficient, due to the high shaft flexibility. 

A new rotor was analyzed in which the shaft stiffness between 
sections 4 and 5 was increased from 0.62 to 1.0 inch. This stiffen- 
ing of the shaft increased the strain energy distribution from 27% 
to 60%. An analysis wa:> also performed with the ^‘jl 1 ey mass attached 
to the shaft. There will only be one critical speed present even 
with the pulley mass included. It was also determined that the 
pulley mass and self-aligning ball bearing will have little effect 
on the dynamic characteristics of the eddy-current test rotor. 

From an examination of the rotor mode shape, it is seen that 
the rotor system may be approximated by a single mass Jeffcott 
rotor using the rotor modal characteristics. The rotor may be 
balanced by single plane corrections at section B, the Inertia mass 
location. The damper may be therefore examined under a large range 
of unbalances and rotor eccentricities to examine nonlinear effects 
In the eddy-current damper. 

In the determination of the damper characteristics, it Is de- 
sirable to have noncontacting Inductance probes to monitor the shaft 
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motion and also strain gauges mounted on the retainer squirrel cage 
spring to determine forces transmitted to the eddy-current damper. 
By means of the shaft displacements end direct determination of 
forces transmitted to the eady-currert damper, the damping and 
stiffness characteristics of the system may be oetermined for a 
range of speeds and unbalances. 
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VIII. DISCUSSION AND CONCLUSIONS 


The object of this investigation was to study the feasibility and 
characteristics of a passive eddy-current damper for application in a 
cryogenic pump. In the design of the cryogenic pump with duplex ball 
bearings, it is extremely difficult to incorporate sufficient damping 
into the system to control the critical speeds or the occurence of self 
excited whirl instability. Flexible supports with coolant friction 
damping have not proven to be successful in cryogenic pumps. 

The concept of the passive eddy-current damper is that the outer 
nonrotating brace of the rolling element b~ supported by a nonferrous 
disk which moves in a permanent magnet field. The eddy-current damper 
has been used for commercial applications such as meters and instruments 
in which a small amount of damping is required. 

The theoretical predictions for the amount of damping generated 
is proportional to the power loss generated by the induced eddy-currents 
created in the damper material. In a simplified eddy-current damper 
analysis, the damping coefficient is proportional to the magnetic flux 
density squared, the damper volume and conductivity. With the new rare 
earth magnetic materials available, a substantially larger field can be 
achieved than with the conventional Alinco magnets. The operation at 
cryogenic temperatures should cause a substantial improvement over 
ambient conditions due to the increase in conductivity of the material. 
There is possibly an off setting factor in that the full damping 
effectiveness is not achieved due to the skin penetration effect of AC 
currents, the higher the frequency of oscillation, the less the penetration 
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depth of the eddy-currents into the conductor. This means that the full 
volumn of the conductor is not available as a damper material. This 
effect needs more experimental testing to be verified. 

Preliminary tests conducted at NASA and at the University of Virginia 
showed a substantial increase in damping characteristics under cryogenic 
temperatures* However the elementary tests conducted at the University 
of Virginia did not simulate the high frequency operation that would be 
encountered with a cryogenic pump operating in a 30,000 rpm speed range. 

The characteristics of an eddy-current damper were analysed from 
several elementary standpoints and also by means of a finite element 
two dimensional analysis program using the electric potential vector 
approach. By means of this finite element analysis, a very complex 
geometry may be computed. However one difficulty in the theoretical 
calculation of the eddy-current damping is the realistic determination 
of the effective magnetic flux density, effective material volume and 
possibly the effect of the depth of penetration. An experimental program 
is necessary to arrive at empherical relationships. Because of these 
problem areas, it was difficult to assign an exact computation of the 
damping generated by the finite element program. 

In general it is felt that the concept of the passive eddy-current 
damper is feasible in a cryogenic pump and further experimental testing 
on this system should be conducted on a high speed simulation model. 

The passive eddy-current damper has the distinct advantage in that it 
has no moving parts and hence is not subjected to fatigue or wear. It 
may also be possible to incorporate into the eddy-current damper design 
a squeeze film damper to generate initial damping characteristics based 
on shear of the fluid film. 
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APPENDIX II 


Finite Element Computer Program and Sample Eddy-Current Problem 

A computer program, Eddy 2, was used to calculate the total power 
loss in an eddy-current damper model. The model was divided into triangular 
elements (see Figure 1) with the power loss from each element summed to 
find the total power loss. 

The total number of triangular elements (Nl) must be input along with 
the total number of nodes, or corners of the triangular elements (N2). The 
number of zero boundary points (N3) is the number of nodes in which the 
elements adjacent to that node have a magnetic field intensity (HO) equal 
to zero. Each triangular element should consist of only one material, 
although a given type of material may exist in many different elements 
throughout the model. N4 is the total number of different materials in 
the model . 

It is also necessary to input the nodes, or vertices [IA(I), JA(I), 
and MA(I)] and the material number [LA(I)] for each of the Nl elements. 

Using some reference point as the origin, the X and Y coordinates of each 
of the N2 nodes, or vertices, must also be input. Additionally, for each 
of the N3 zero boundary points, input the node or vertex number and the 
magnetic field intensity at that node. 

The final data input is for the material properties. For each of the 
N4 materials, the relative permeability of the material (XIRON), the 
magnetic field intensity (HO), and the conductivity (SIGMA) must be added 


A2-1 








to the data input. Note that the magnetic field intensity (HO) is a 
complex number. For the calculation of the current density (CUR) in the 
program, the material relative permeability and conductivity has been set 
to 1. 

Additionally, the frequency is read into the program within the 
program itself, and is currently set equal 60 Hz. 

In Figure 1, the model is divided into 72 triangular elements with 
49 nodes. Three different sets of material properties are present, witn 
the magnetic field intensity equal to zero in all elements closest to the 
boundary of the model except elements 30 and 43. Therefore, there are 24 
nodes which are not in contact with an element which has a magnetic field 
intensity greater than zero. 

The elements were designated by using the vertices and material type. 
Using Node 1 as the origin, the normalized X and Y coordinates of each 
vertex are given using the axis of symmetry as the X-av<s. The 24 zero 
boundary points were all given a magnetic field intensity value of zero. 
Finally, the material properties are read in using a relative permea- 
bility equal to 1.0 and a conductivity equal to 1.0 in all three cases. 

The magnetic field intensity increases from 0 + Oi in the outer elements, 
to 6 + Oi in the intermediate elements, and to 10 + Oi in 8 of the 
central elements. These values are all entered using the MRS system of 
units. 

The listing for the program Eddy 2 along with the printout of the 
input and output data, which shows the total power loss, is given next. 
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PROGRAM EDDY 2 

FIRST ORDER FINITE ELEMENT SOLUTION 
FOR ELECTROMAGNETIC FIELD PROBLEMS 
CF X , Y PLANAR GEOMETRIES INCLUDING 
EDDY CURRENTS IN CONDUCTING PARTS, 
BUT EXCLUDING THE SOURCE REGION. 

EVP METHOD 
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THIS PROGRAM USES ONLY TRIANGULAR FINITE ELEMENTS 

ALL DIMENSIONS AND MATERIAL SPECS 
ARE IN RATIONALIZED M.K.S. UNITS 

COMPLEX S, P, T, CUR, HO 

DIMENSION S < 1 00 , 1 00 > , P < 1 00 > , T < 1 88 ) , CUR < i uA > , HOC 108) , DELTAMC 106;* 
COMPLEX SI I , SI J, SIM, SJJ, SJM, SMM, TT , T I , T J, TM, FACTOR 
DIMENSION IRC 200 ) , JR <208) , MAC 200) , LA <200) , XNUC200) , PE 1O0> 
DIMENSION XU00) , Y< 108) , IXC58), XIRONC 18) , SIGMAC 8> 

CHARACTER*80 HEAD1 , HEAD2 
POWER-8. 8 
AI-BB-BXX-BYY-0. 0 
INITIALIZE ARRAYS 
DO 990 1=1,109 
DO 998 J*l, 100 
SCI, J>*<0. 0,0.0) 

998 CONTINUE 
990 CONTINUE 

DO 997 1*1, 108 
PU>*<0. 8,0.0) 

TU)«(0. 0,0.0) 

CURCI>*<0. 0,0.0) 

997 CONTINUE 

DO 996 1*1 , 200 
I A ( I )*0 
JAC I )=0 
MAC I ) =0 
LA C I ) *0 
XNU < I ) *0 . 0 
996 CONTINUE 

DO 995 1*1 , 100 
X < I ) =0 . 0 
PE C I )*0. 0 
DELTAMC I )*0. 0 
HOC I )*C0. 0, 0. 0) 

YCI)*0.0 
995 CONTINUE 

DO 994 1*1,50 
IXC I )*0 
S94 CONTINUE 

DO 993 1*1, 10 
XIRONC I )*0. 0 
SIGMAC I >*0. 0 
993 CONTINUE 

NOPT *0 

READ MODEL CONTROL DATA 

IB IS THE FILE LINE NO. 

Nl-TOTAL NO. OF ELEMENTS 
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N2“T0TftL NO. OF NODES 
N3=N0. OF ZERO BOUNDfi.'Y PTS. 
N4=N0. OF MATERIALS 
READ (5, 778>HEAD1 
REALc.5, 778>HEAD2 
778 FORMAT ( A80) 
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READ<5,*)IB,N1,N2,N3, N4 

ELEMENT CONNECTIVITY 
AND MATERIAL ID NO. 


IA, JA,MA*VERTICES OF TRIANGULAR ELEMENTS 
LA=MATERI AL ID NUMBER 
DO 12 1*1, N1 

12 READ(5, *> IB, I ft C I > , JA< I > , MA< I > , LA Cl ) 

READ NODE POINTS 

X, Y=COORBI NATES OF VERTEX 
DO 13 I = 1 , N2 

13 REABt5,*>IB,X<I>, Y<I> 

BOUNDARY PTS. & BOUNDARY POTENTIALS 
I X=N0DAL NUMBER OF ZERO BOUNDARY POINTS 

DO 14 I * 1 f N3 

READ <5, *> IB, IX<I>,TEMP 

14 P < I X < I ) > =CMPLX < TEMP , 0 . ) 

READ MATERIAL PROPERTIES 
XIRGN=RELATI VE PERMEABILITY 
HO=MAGNETIC FIELD INTENSITY 
CUR=CURRENT DENSITY 
SIGMR=CONBUCTIVITY 
FREQ*FREQUENCY OF EXCITATION IN hZ 
OMEGA*ANGULAR FREQUENCY 

XMU0*4. *3. 14159*d.0E-7> 

FREQ*60. 

0MEGA*2. *3. 141 59*FREQ 


DO 120 1=1, N4 

READ (5, *> IB, XIRONt I > , H0< I > , SIGMA< I > 

CUR( I >=-CMPLX<8, 1 > ?*QMEGA*S I GMA < I >#H0< I > 

IFCXIRONC I ) . EQ. 0)XIRONC I >=1 . 

IF(.SIGMA< I ) . EQ. 0)SIGMA< I > = 1 . 0E-6 
120 XIR0N< I) = l./<XIRONa^XMUO> 

WRI TE<6, 778 ) HERD 1 
WRITERS, 778 )HEAD2 
WRI TE<6, 201 > 

201 FORMAT < 1H ,//,lX, H EBDYfl IX, "ELECTROMAGNETIC FIELD SOLUTION 1 ', 
1/,1X," OF X , Y PLRNAR GEOMETRIES" , IX, " BY THE FINITE ELEMENT 
2METH0D" ) 

WRI TE < 6 , 102 )N1 

102 FORMATC1H , 2X , 11 TOTAL NO. OF ELEMENTS *",I4> 

WRITEC6, 1 03 > N2 

103 FORMAT < 1H ,2X," TOTAL NO. OF NODES »“,I4> 

WRITEC6, 1 04 >N3 

104 FORMAT C 1 H ,2X," NO. OF ZERO BOUNDARY PTS *",I4> 

WRITE <6, 1 05 > N4 

105 FORMATOH ,2X," NO. OF MATERIALS *",I4, 

1 //, *' ******************************************* ****** " ) 

WRITERS, 999 > 
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999 FORMAT < 1H ,////////, "MATERIAL PERMEAB ILI TV , CONDUCTIVITY «. SOURCE 
1 CURRENT DENSITY",//) 

DO 12:00 1*1, N4 
IB* I 

WRITE (6, 1201 ) IB, XI RON < I ) , SIGMA< I ) , CUR< I ) 

1201 FORMAT <15, 4E 15.4,//) 

1200 CONTINUE 

CONSTRUCT MATRIX 

DO 5 1*1, N1 
xnu<I)=xironi:la<I)> 

FIND GEOMETRIC COEFFICIENTS 

BI*Y<JA<I>>-Y<MA<I)) 

BJ=Y<MA< 1 > >-Y< I A< I'/j 
BM*YtIAU)>-Y<JA<I>) 

C!=X(MA<I))-X<.JA(I)) 

CJ*X<IA',I))-X(MA<I)> 

CM*X<JA(I))-X<IA<I)) 

D 1 *X < JA < I > ) * Y v MA < I i :« 

D2*X<MA<I))*Y<JA<I)'- 
D3*X< 1A< I ) )*BI 
D4=Y< IA< I ) )*CI 

COMPUTE AREA OF TRIANGLE 

DELTA=A S<<Dl-D2+D3+D4>/2. ) 

DELTAM< I )=DELTA 

FACTOR*CMPLX<0. , 1. )*OMEGA*SIGMA<LAi' I ) >*DELTA/12. 

COMPUTE MATRIX ELEMENTS 

BB=<BI*BI+CI*CI )/<4. *DELTR; 

BC=<BI*BJ+CI*CJ)/<4.*DELTA> 

BD»<BI*BM+CI*CM)/<4.* DELTA) 

CC=<BJ*BJ+CJ*CJ)/(4.*DELTA) 

CD=<BJ*BM+CJ*CM)/<4.*DELTA) 

DD=<BM*BM+CM*CM)/<4.*DELTA) 

SI 1*2. *FACTOR+BB*XNU<' 1 ) 

SIJ=FACTOR+BC*XNU<I) 

SIM*FACTOR+BD*XNUa ) 

SJJ*2. *FACTGR+CC*XNU( I ) 

SJM*FACTOR+CD*XNU( I ) 

SMM*2. *FACTOR+DD*XNU( I ) 

ASSEMBLE MATRIX ELEMENTS 

S<IA<I>, IA<I))«S<IAU>, IA(1))+SII 
S<IA<I), JA<I))*S<IA<I), JA<I))+SIJ 
S<IA<I),MA<I))*S<IA<I),MA<I))+SIM 
S<JA<I),JA<I))*S<JA'.I),JA(I))+SJJ 
S<JA<I>,MA<I))*S<JAU. ,MAi:i))+SJM 
S<MA<I),MA<I) >*S<MA*. I . ,MA<I J)+SMM 

S< JA< I ) , IA<I))*S<IA<I), JA<I>) 

S(MA< I ) , IA<I))*S<IA<I),MA<I)) 

S<MA<I),JA<I))*S<JA<I),MA(I)) 

COMPUTE FORCING FUNCTION FOR CURRENT REG t ON 

TT*CUR<LA< I ))* DELTA/3. 

C 
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ASSEMBLE FORCING FUNCTION FOR CURRENT REGION 


T I =TT 
T J*TT 
TM*TT 

Tafl(i))*Taft<i)) i tt 
I tJflU)>*T<JfltI)>*TT 
T<MA<iI))*T<MA<I>) + T7 

PRINT ALL DATA 

I FCNOPT. EQ. 0>GQ TO 5 
WRITE(6, 777) I 

777 FORMAT < 1H , / , “ELEMENT NUMBER: " , 1 2 > 

WRI TE< 6 , 1000)BI,BJ,BM,C1,CJ,CM 
1600 FORMAT < 1H ,/',lX," GEOMETRIC COEFFICIENTS" , 5X, "BI “ , 10X, 

1 "BJ“ , 10X, "BM", 10X, "Cl", 10X, "CJ", 10X, "CM" , /, 6E12. 4) 

WRITE<6, 100DDELTA 

1001 FORMAT < 1 H IX," AREA OF TRIANGLE*" , E12. 5) 

WRITERS, 1002)BB,BC,BD,CC,CD,DD 

1002 FORMATOH ,/v,lX,“ MATRIX ELEMENTS" 5X , " BB ", 1 0X , " BC" , 10X, " ED " , 
U0X, "CC", 10X, "CD", 10X, "DD" , ' , 6E12. 4) 

WRITEC6, 1 003 > 

1003 FORMAT < 1 H INDIVIDUAL ELEMENT MATRIX".) 

WRITE <6, 1004>S:i,SIJ,SIM 

WRITE <6, 1004>SJ J,SJJ,SJM 
WRITE<6, 1004)SIM,SJM, SMM 

1004 F0RMAT<1X,6E12.4> 

WRITE<6, 1005) 

1005 FORMAT U H IX," FORCING FUNCTION") 

WRITE<6, 1 008 >T I , T J, TM 

1008 F0RMATdX,6F15.4) 

WRITEC6, 1006) 

1006 FORMAT <.//////} 

5 CONTINUE 

WRITE<6, 1019) 

1019 FORMAT < 1H ,/v,lX," FULL FORCING FUNCTION") 

DO 1016 K= 1 , N2 
1016 WRITE<6, 300? )K, TOO 
3007 F0RMAT(I5,E12.4,E12.4) 

1050 CONTINUE 

ASSEMBLE FORCING FUNCTION TO INCLUDE 
APPLIED POTENTIALS - IF ANY 

DO 151 IC=1,N3 
DO 151 J*1,N2 

151 T<J)"T<J>-<S<J,IX<IC)>*P<IX<IC>>> 

MODIFY FORCING FUNCTION ACCORDING TO 
BOUNDARY CONDITIONS 

DO 2000 IC«1,N3 
DO 2000 1*1, N2 
IF<I- IX Cl 0)2000, 2003, 2000 
2003 T<I)*P<IX<IC>> 

MODIFY [S3 MATRIX ACCORDING TO 
BOUNDARY CONDITIONS 

DO 2002 >1,N2 
S< J, I )*0. 

IFCJ.EQ. I>SU, J)*l. 

S< I , J)*S< J, I ) 
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2O02 CONTINUE 

SET FORCING FUNCTION ROW = 
APPLIED POTENTIAL - IF ANY 
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IFCI.EQ. IXCIC) >T<rn=P( IXt TO) 

2000 CONTINUE 


MORE PRINT 
WRITEC6, 1011) 

1011 FORMAT < 1 H ,11V 1 MODIFIED FORCING FUNCTION" ) 
DO 1012 K= 1 , N2 

1012 WRITE<6,300?>K, T<K) 


1009 CONTINUE 

CALL SOLVE ( P , T , S , N2 ) 

WRITE POTENTIALS 
WRITERS, 111) 

111 FORMAT OH , IX, “ELEMENT 1 2X , " COMPLEX FLUX VECTORS", 
+ /,2X,78<"*">> 

DO 16 I = 1 , N2 
16 WRITERS, 1?) I,PU* 

1? FORMATOH , 5X, 12, 9X, " < " , El 5. 6, " , - , E15. 6, " > " > 


EVALUATE VECTOR POTENTIALS 

WRITEC6, 1020) 

1020 FORMAT < 1H , 1 7X, "VECTOR POTENTIALS") 

1021 FORMATOH , 6X, " EXX " , 1 3X, " BYY " , 1 3X , "BED" ) 

WRI TE<6, 1022) 

1022 FORMAT < 1H , " *** ****** ********************************************** 

+****************" > 

WR I TE ( 6 , 222 > 

222 FORMAT < 1H , IX, "ELEMENT #" , 7X, " JXM" , 12X, " JYM" , 12X, " JM" , 

+8X, "LOCAL POWER") 

C 

DO 8001 I = 1 , N 1 

AI*CX(JA<I))*Y<MA<I>)“X(MA<I)>*Y<JA<I))> 

BI=<Y<JA<I))-Y<MA<I>)> 

CI*<X<MA<I))-X<JA<I>>) 

CJ*XUA<I))-X<MR<I>) 

CM*X< JAC I))-X<Ifl<I)> 

BJ*Y<MA<I>)-Y<IA<I>) 

BM*Y<IA<I>)-Y<JAU)> 

C 

BYY«REflLCBI#Pafi<I>>+BJ*PUA<n> + BM*P<Mfl<I>)> 

BYY*-BYY/<2.*DELTAM<I )) 

BXX=REAL<CI*P<IA<I ))+CJ*POfl< I >)+CM*P(MA< I ) ) ) 

BXX**BXX/<2. * DEL TAM (I ) ) 

BB«<BXX**2+BYY**2> 

PE< I >*BB*DELTRM< I >/SIGMA(LA< I ) ) 

WRITERS, 333) I , BXX, BYY , BE, P£<. I > 

333 FORMAT < 1 H , 4X, 12, 4X, 4<E15. 6) ) 

POWER*POWER+PE< I ' 

8001 CONTINUE 

WRI TE <6, 888 ) POWER 

888 FORMAT < 1 H 49X , 29 47X, " TOTAL POWER * " ,E1 5 . 6 ) 

C 

STOP 

END 


i I 


\ 

i 



o o 



c 

SUBROUTINE SOLVE<P, T, S, N2> 

C 

C SOLVE SIMULTANEOUS EGfUA T IONS 

BY GAUSSIAN ELIMINATION METHOD 

COMPLEX P,T,S 

DIMENSION Pa00),T<100),S<:00, 100) 
COMPLEX FF,FFX 
FF=(3. 0,0.0) 

FFX=<0. 0,0.0) 

K = 0 

i k'X=0 

!\ ' c 

DO 111 1=1, N2 
P< I ) = T< I ) 

111 CONTINUE 
C 

M*2 

NN=N2-1 

C 

DO 1 1=1, NN 

C 

DO 2 J*M,N2 

. • FF=-S(J, I)/S(I, I) 

* P<J)=P< J)+FF*P< I ) 

. : C 

Li DO 3 K=M,N2 

- -J ? S<J,K>=S<J,K)+FF*SU.K> 

* ! C 

! 2 CONTINUE 

* 1 M=M+ 1 

C 

P<N2^=PCN2)^S<N2,N2) 

M=N2- 1 

, C 

* ■ DO 100 1=2, N2 

| J=N2- I + 1 

* ! FFX*P<J) 

' ; NN*N2-1 

i C 

DO 200 K=M,NN 
. j KX*K+1 

' 200 FFX=FFX-S< J,KX)*P<KX) 

* ; C 

* j M*M- 1 

j 100 P< J)=FFX^S< J, J) 

■ ! C 

, , RETURN 

S j END 

I;, j END OF FILE 
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2 POLE COBALT - RARE EAPTH MAGNET 
APPROX. 1 IN. DIA. / .3 75 IN. HIGH, ASSUME 0.25 IN. 
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EDDY** 

ELEC TROMr CHET IC FIELD SULUT ION 
OF FLAHftR GEOilETRIES 

BY THE FINITE ELEMENT METHOD 

TOTfiL NO. OF ELEMENTS = 72 

TOTAL NO. OF NODES = 49 

NO. OF ZERO BOUNDARY PTS = 24 

NO. OF MATERIAL 3 - 3 


\ ; MATERIAL F ERMEAB I L IT Y , CONDUCT I V I T Y & SOURCE CURRENT DENSITY 

k t 

; * 1 . 7958E+06 . 1000E + 01 0. 0. 


2 . 7958E+0E . 1 006E+0 1 0. 


2Z62E + 04 


.7958E+0* 


. 1O00E+01 O. 


- . 37701: +04 


I 
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FULL 

FGFCING 

FUNCTION 
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0 . 
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0. 
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0. 

13 

0. 

-.2011E+04 

14 

0. 

2953E+04 

15 

0. 

1131 E+04 

16 

0. 

0 » 

17 

0. 

5027E+04 

18 

0. 

-.2765E+04 

19 

0. 

~ . 754OE+03 

20 

0. 

0. 

21 

0. 

2727E-r04 

22 

0. 

235CE+04 

23 

0. 

Ill 2E+04 

24 

0 . 

0. 

25 

0 • 

0. 

26 

0. 

- . 754OE+03 

27 

0. 

1 093E+04 

23 

0 # 

d. 

29 

0 * 

8 . 

30 

0. 

0. 

31 

0. 

0. 

32 

0. 

1 093E+04 

33 

0. 

0. 

34 

0. 

0. 

35 

0. 

Ill 2E+04 

36 

0, 

2350E+04 

37 

0. 

0. 

38 

0. 

7540E+O3 

39 

0. 

2765E+04 

40 

0. 

0. 

41 

0. 

- . 1 131E+04 

42 

0. 

-.2953E+04 

43 

0. 

0 » 




1 


i 

A 


;S 

i 





/: 

/I 

M 





.1 

> 

* r 

-j ' 


i 

, t 


* I 
? 


44 

0. 

1 508E+04 

45 

0. 

0. 

46 

A 

-. 1 500E+84 

4? 

0. 

0. 

48 

6. 

0 , 

49 

'!. 

0. 


ELEMENT # COMFlEX FLUX VECTORS 


ORIGINS FW£ ® 
OF POOR QUALITY 


* * * * * i 




*** *-►•* 


»£. 


o. , 

u. > 

3 

f 

0. , 

0. > 

4 

i 

0. S 

0 . :■ 

5 

< 

2 3 S 6 9 5 E - 0 5 , 

-.21 8690E-02 > 

6 

i 

0. 

0. ; 

*7 

1 

< 

479827E-05, 

- . 4671 06E-02 ) 

8 

< 

327412E-05, 

- . 3 G 7 0 1 8 E - 6 2 ’> 

9 


0. , 

0 . > 

10 


- . 78 3586E-05 , 

-. 734258E-02; 

11 

s 

-.431638E-65, 

- . 442330E-62 ) 

12 

K 

0. , 

0 . > 

13 

< 

-.95922OE-05, 

-.891 762E-82) 

14 


- . S34467E-05 , 

- . 76936 1 E-02 > 

15 

< 

470924E-05, 

-.411 365E-G2 > 

16 


0. , 

0. > 

17 

( 

Ill 269E-04 , 

-. 104141E-01) 

18 


-.91S324E-05, 

-. S43161E-02.) 

19 

< 

- . 505058E-85 , 

- . 4363 32E-62 ) 

20 

( 

0. , 

0. > 

21 

f" 

106836E-04, 

-. 955S86E-02) 

22 

< 

- . 865691 E-05 , 

- . 7781 35E-02 > 

23 

< 

- • 352464E-05 , 

3039O6E-02) 

24 

( 

0. , 

0. > 

25 

( 

0. , 

0. / 

26 

< 

-. 101566E-04, 

883125E-02.) 

27 

< 

- . 82O247E-05, 

710439E-©2> 

-. 

< 

0. , 

0. > 

29 


0 « 

0. > 

30 

< 

0. , 

0. 

31 

< 

0- 

0. > 

32 

C 

- . S20247E-05, 

710439E-02> 

33 

( 

0. , 

0. > 

34 

c 

6. , 

0. > 

35 

< 

- . 352464E-05, 

-. 303906E-02) 

36 

< 

- • 865691 E-05, 

- . 7781 35E-02 > 

37 

< 

0. , 

0. > 

38 

< 

- . 505058E-05, 

-.436332E-02) 

39 

< 

-.91 6324E-05 , 

- . 843161E-02 > 

40 

( 

0. , 

0. > 

41 

< 

47G924E-05, 

-. 41 1365E-02> 

42 

( 

- . 83446 7E-05, 

-. 76936 IE-02) 

43 


8. , 

0. > 

44 

< 

■* . 43 1638E-05, 

- . 442330E-02 > 

45 


0- , 

0. ) 

46 

< 

- . 32741 2E-05 , 

- . 3070 1 8E-02 ) 

47 


G. , 

0. > 

48 

< 

0. , 

G. > 

49 

i 

0 . , 

0. / 


VECTOR POTEHTIRLS 


**********- 

**************** 

**************** 

***************************** 

ELEMENT # 

JXM 

JVM 

Jh 

LOCAL POWER 

1 

0. 

0. 

0. 

0. 

2 

0. 

0. 

0. 

0. 

3 

0. 

23869 5 E-05 

. 5b'3753E-l 1 

. 404525E- 1 1 

4 

1 29O24E-05 

- . 554305E-06 

. 197254E-1 1 

. L 82460E- 1 1 

5 

1 229 1 6E-G5 

-.2411 32E-05 

« 732529E-1 1 

.454168E-11 

6 

1 29824E-05 

- . 24S707E-05 

. 785024E- 1 1 

. ' 7 26 1 47E- 1 1 

7 

314819E-65 

- . 13 537 2 E-05 

. 1 1 7437E- 1 0 

. I51-.S71E-1 1 



A2* 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


■ 1 

8 

122 916E-85 

- . 303759E-05 

. 1 07378E- 1 0 

.i;65741E-11 

9 

1 94632E-65 

- , 2 1 4830E-05 

. 840337E- 1 1 

.•530252E-11 

i 5 

- t 

10 

314319E-05 

- . 246079E-05 

• 1 5 9666E- 1 0 

. 830263E- 1 1 

1 1 

409J93E-05 

1 72073E-05 

. 1 97462E- 1 0 

. 1 18477E-10 

1 V 

12 

-, 1 24 753E-85 

- • 1 75634E-05 

. 464 1 C7E- 1 1 

. 2328S3E-1 1 

; !_ 

13 

-. 1 94632E-05 

-.24551 3 E- 05 

. 981583E-1 1 

.* ? 36187E-11 

* 1 

fc 

14 

-.363543E-05 

1 61 058E-05 

. 1 53 1 03E- 1 0 

. ' 7 965 1 6E- 1 1 

. V 

15 

-.409698E-05 

184135E-05 

• 20 1 758E- 1 0 

. 1 21055E- 1 0 

1* * 

16 

-. 416747E-85 

179201E-05 

• 26579 1 E- 1 0 

. L 1 6272E- 1 0 

f 

17 

-. 1 24753E-05 

1 53472E-05 

. 391 170E-1 1 

. L95585E-1 1 

■ 

IS 

-. 1 96368E-05 

-.818573 E-06 

. 4526 1 0E- 1 1 

. 226385E- 1 1 

t 

19 

363543E-05 

-.818573 E-06 

. 1 33S64E- 1 0 

• 694321 E-l 1 

f ' 

t 

20 

-.411 266E-05 

341344E-06 

. 1 73385E- 1 8 

. 351525E-1 1 

\ » 

21 

-.41 6747E-05 

341344E-06 

. 1 74843E- 1 0 

. 9S7863E- 1 1 

J * ' 

22 

-.321 593E-35 

141 545E-05 

. 12352 IE- 10 

. 969644E- 1 1 

L • 

23 

202605E-&5 

• 44396 1 E-06 

. 430197E-1 1 

. 2 1 5099E- 1 1 

“ y 

24 

1 96368E-05 

. 5O6330E-06 

. 4 1 1241E- 1 1 

. 205629E- 1 1 

J 

25 

-.51 3228E-85 

. 506330E-66 

.265966E-10 

. 1 32983E- 10 

C 

26 

-.41 12S6E-05 

. 1 52595E-05 

. 192425E- 10 

• 962 1 25E- 1 1 


27 

321593E-65 

. 152595E-05 

. 1 26772E- 1 0 

. 995 1 57E- 1 1 

j.' 

28 

352464E-05 

. 200904E-05 

. 1 64593E- 1 6 

. 322966E-1 1 


29 

0. 

0. 

0. 

0. 


30 

1 98048E-85 

. 526331E-06 

. 4 1 9932E- 1 1 

. 209966E- 1 1 

• « 

31 

-. 202605E-05 

. 478364E-06 

. 433770E- 1 1 

.205851E-11 

« * a 

32 

5 13228E-05 

. 478364E-86 

. 265691E-10 

. 126203E-10 

| 

33 

914237E-05 

. 469952E-05 

. 1 05668E-09 

. 396256E- 1 8 

► ! 

34 

1 40985E-05 

. 46 9952E-05 

. 240731E- 10 

. 902743E- 1 1 

>' *; 

35 

- • 352464E-05 

0. 

. 1 2423 1 E- 1 0 

. 279519E-1 1 

fc 

i 

36 

0. 

0. 

0. 

0. 

* | 

37 

0. 

0. 

0. 

0. 

t 

38 

• 352464E-05 

0. 

. 1 24231 E- 1 0 

• 279519E-1 1 

5" ~ ! 

39 

. 1409S5E-05 

. 469952E-05 

. 24073 IE-10 

• 902743E- 1 1 


40 

, 9 1 42 37E-05 

. 469952E-05 

. 1 05668E-09 

.396256E-10 

: 

41 

. 513228E-65 

. 478364E-06 

. 26569 1 E- 1 6 

. L 26203E- 1 0 


42 

. 202605E-05 

. 478364E-06 

„ 433370E- 1 1 

• 205851E-1 1 


43 

. 198048E-05 

. 526331 E-06 

. 4 1 9932E-1 1 

. 209966E- 1 1 

. 

» 

44 

0. 

0. 

6. 

0. 

■* 

45 

352464E-05 

• 200904E-05 

. 164593E-10 

. 322966E- 1 1 


46 

\ 321 693E-05 

. 1 52595E-05 

. 1 26772E- 1 0 

. 995 1 57E- 1 1 

j • 

47 

.411 266E-05 

. 1 52595E-05 

. 1 92425E- 1 0 

. 962125E-1 1 


48 

• 5 1 322SE-05 

. 506330E-06 

. 265966E- 1 0 

. L 32983E- 1 0 


49 

. 1 96368E-05 

. 506330E-06 

. 41 1241E- 1 1 

• 205620E- 1 1 

' t 

58 

. 202505E-05 

. 443961E-06 

. 430 1 97E- 1 1 

. 21 5099E- 1 1 

• t 

i 

51 

.321 693E-05 

141 545E-05 

. 1 2 352 1 E- 1 0 

. 969644E- 1 1 

I 

1 s 

52 

. 416747E-65 

341344^-06 

. 1 74843E- 1 0 

. 987S63E- 1 1 

\ 

53 

.411 266E-05 

-. 341344E-06 

. 1 70305E- 1 0 

. 85 1 525E- 1 1 

t 

54 

. 363543E-C5 

-.81 S573E-06 

. 1 3S864E- 1 9 

. 69432 1 E- 1 1 

i 

55 

. 196368E-05 

81S573E-06 

. 4526 1 0E- 1 1 

• 226305E-1 1 


56 

. 1 24753E-05 

-. 1 53472E-05 

.3911 70E- 1 1 

. 1 95585E- 1 1 


57 

. 4 16747E-05 

-. 179201E-85 

. 20579 i.E- 1 0 

. 1 16272E-10 

r 

58 

. 409698E-05 

1841 35E-05 

. 20 1 758E- 1 0 

. L 2 1 055E- 1 0 

! 

59 

. 363543E-85 

-. 161 058E-O5 

. 1581 03E-1 0 

. ' 7 905 1 6E- 1 1 


60 

. 1 94632E-05 

- • 24551 3E-05 

• 981583E- 1 1 

• ? 361 87E-1 1 


6 1 

. 1 24753E-85 

1 75634E-05 

. 464 1 07E- 1 1 

. 232053E- 1 1 


62 

• 409S98E-05 

-. 1 72073E-05 

. 1 97462E- 1 0 

. 1 1 8477E- 1 0 


, 63 

. 3 1 4 8 1 9 E - 8 5 

- • 246079E-05 

. 1 59666E- 1 0 

. 330263E- 1 1 


: 64 

• 1 94632E-05 

-. 214830E-U5 

. 840337E- 1 1 

. 630252E- 1 1 


65 

. 122916E-05 

- . 363759E-05 

. 1 0737SE- 1 0 

. 66574 1 E- 1 1 


! 66 

. 314819E-05 

135372E-05 

. 117437E-10 

. 6 1 067 1 E- 1 1 


i 67 

. 1 2291 6E-05 

-.2411 32E-05 

• 732529E- 1 1 

• 454 1 68E- 1 1 


68 

. 1 29024E-G5 

-.2 48707E-05 

. 785024E- 1 1 

. ? 26 1 47E- 1 1 

r 

69 

. 1 29024E-85 

- . 554805E-06 

. 197254E-11 

. 1 82460E- 1 1 


70 

0. 

238695E-05 

. 569753E- 1 1 

• 404525E- 1 1 

nr - 

71 

0. 

O. 

0. 

0. 

r 

72 

0. 

0. 

0. 

0. 

a 
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. 5O1995E-09 


I® 



